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A new approach to the single-band Hubbard model is described in the general context of many- 

fl^ . body theories. It is based on enforcing conservation laws, the Pauli principle and a number of 

crucial sum-rules. More specifically, spin and charge susceptibilities are expressed, in a conserving 

approximation, as a function of two irreducible vertices whose values are found by imposing the 

. ■ local Pauli principle (j^) — (n|) as well as the local-moment sum- rule and consistency with the 

equations of motion in a local- field approximation. The Mermin- Wagner theorem in two dimensions 

is automatically satisfied. The effect of collective modes on single-particle properties is then obtained 

by a paramagnon-like formula that is consistent with the two-particle properties in the sense that the 

1-^ I potential energy obtained from TrEG is identical to that obtained using the fluctuation-dissipation 

(~| , theorem for susceptibilities. Since there is no Migdal theorem controlling the effect of spin and 

Q ' charge fluctuations on the self-energy, the required vertex corrections are included. It is shown 

O , that the theory is in quantitative agreement with Monte Carlo simulations for both single-particle 

and two-particle properties. The theory predicts a magnetic phase diagram where magnetic order 

persists away from half-filling but where ferromagnetism is completely suppressed. Both quantum- 

^ , critical and renormalized-classical behavior can occur in certain parameter ranges. It is shown that 

QQ ■ in the renormalized classical regime, spin fluctuations lead to precursors of antiferromagnetic bands 

OO \ (shadow bands) and to the destruction of the Fermi-liquid quasiparticles in a wide temperature range 

above the zero-temperature phase transition. The upper critical dimension for this phenomenon is 

04 ' three. The analogous phenomenon of pairing pseudogap can occur in the attractive model in two 

^^ \ dimensions when the pairing fluctuations become critical. Simple analytical expressions for the self- 

y_^ ■ energy are derived in both the magnetic and pairing pseudogap regimes. Other approaches, such 

as paramagnon, self-consistent fluctuation exchange approximation (FLEX), and pseudo-potential 

parquet approaches are critically compared. In particular, it is argued that the failure of the FLEX 

a approximation to reproduce the pseudogap and the precursors AFM bands in the weak coupling 

regime and the Hubbard bands in the strong coupling regime is due to inconsistent treatment of 
vertex corrections in the expression for the self-energy. Treating the spin fluctuations as if there 
pH ' was a Migdal's theorem can lead not only to quantitatively wrong results but also to qualitatively 

Q I wrong predictions, in particular with regard to the the single-particle pseudogap. 
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C^ ; I. INTRODUCTION 

Understanding all the consequences of the interplay between band structure effects and electron-electron interactions 
remains one of the present-day goals of theoretical solid-state Physics. One of the simplest model that contains 
the essence of this problem is the Hubbard model. In the more than thirty years [|l| since this model was 
formulated, much progress has been accomplished. In one dimension, [^ IJ] various techniques such as diagrammatic 
rcsummations, ||^ bosonization, [^ renormalization group Q and conformal approaches |^ |1^ have lead to a 
very detailed understanding of correlation functions, from weak to strong coupling. Similarly, in infinite dimensions a 
dynamical mean-field theory [pTi leads to an essentially exact solution of the model, although many results must be 



obtained by numerically solving self-consistent integral equations. Detailed comparisons with experimental results on 
transition-metal oxides have shown that three-dimensional materials can be well described by the infinite-dimensional 
self-consistent mean-field approach. ||ll| Other methods, such as slave-boson |12| or slave-fermion jlj] approaches, 
have also allowed one to gain insights into the Hubbard model through various mean-field theories corrected for 
fluctuations. In this context however, the mean-field theories are not based on a variational principle. Instead, they 
are generally based on expansions in the inverse of a degeneracy parameter, |14[| such as the number of fermion flavors 
N, where N is taken to be large despite the fact that the physical limit corresponds to a small value of this parameter, 
say N = 2. Hence these theories must be used in conjunction with other approaches to estimate their limits of 
validity. |l3| Expansions around solvable limits have also been explored, flfl] Finally, numerical solutions, |1^ with 
proper account of finite-size effects, can often provide a way to test the range of validity of approximation methods 
independently of experiments on materials that are generally described by much more complicated Hamiltonians. 

Despite all this progress, we are still lacking reliable theoretical methods that work in arbitrary space dimension. 
In two dimensions in particular, it is believed that the Hubbard model may hold the key to understanding normal 
state properties of high-temperature superconductors. But even the simpler goal of understanding the magnetic phase 
diagram of the Hubbard model in two dimensions is a challenge. Traditional mean-field techniques, or even slave- 
boson mean-field approaches, for studying magnetic instabilities of interacting electrons fail in two dimensions. The 
Random Phase Approximation (RPA) for example does not satisfy the Pauli principle, and furthermore it predicts 
finite temperature antiferromagnetic or spin density wave (SDW) transitions while this is forbidden by the Mermin- 
Wagner theorem. Even though one can study universal critical behavior using various forms of renormalization group 
treatments [ISl |19] pffl |21| |22] or through the self-consistent-renormalized approach of Moriya [E3| which all satisfy 



the Mermin- Wagner theorem in two dimensions, cutoff-dependent scales are left undetermined by these approaches. 
This means that the range of interactions or fillings for which a given type of ground-state magnetic order may appear 
is left undetermined. 

Amongst the recently developed theoretical methods for understanding both collective and single-particle properties 
of the Hubbard model, one should note the fluctuation exchange approximation [g^ (FLEX) and the pseudo-potential 
parquet approach. |2^] The first one, FLEX, is based on the idea of conserving approximations proposed by Baym 
and Kadanoff. pO] [27| This approach starts with a set of skeleton diagrams for the Luttinger-Ward functional J28| to 
generate a self-energy that is computed self-consistently. The choice of initial diagrams however is arbitrary and left 
to physical intuition. In the pseudo-potential parquet approach, one parameterizes response functions in all channels, 
and then one iterates crossing-symmetric many-body integral equations. While the latter approach partially satisfies 
the Pauli principle, it violates conservation laws. The opposite is true for FLEX. 

In this paper, we present the formal aspects of a new approach that we have recently developed for the Hubbard 
model [E9[ [ pO| . The approach is based on enforcing sum rules and conservation laws, rather than on diagrammatic 
perturbative methods that are not valid for interaction U larger than hopping t. We first start from a Luttinger-Ward 
functional that is parameterized by two irreducible vertices Ugp and Uch that are local in space-time. This generates 
RPA-like equations for spin and charge fluctuations that are conserving. The local-moment sum rule, local charge sum 

rule, and the constraint imposed by the Pauli principle, (^t ) = ("-t) then allow us to find the vertices as a function of 

double occupancy {n^ni) (see Eqs.(p7|) and (p8|)). Since {n-^ni) is a local quantity it depends very little on the size of 
the system and, in principle, it could be obtained reliably using numerical methods, such as for example Monte Carlo 
simulations. Here, however, we adopt another approach and find (mnj) self-consistently p9| without any input from 



outside the present theory. This is done by using an ansatz En. (40) for the double-occupancy {n^ni) that has been 
inspired by ideas from the local field approach of Singwi et al. ||31|. Once we have the spin and charge fluctuations, the 



next step is to use them to compute a new approximation, Eq.(46), for the single-particle self-energy. This approach 
to the calculation of the effect of collective modes on single-particle properties l^^l is similar in spirit to paramagnon 
theories. 132] Contrary to these approaches however, we do include vertex corrections in such a way that, if T,^^' is our 
new approximation for the self-energy while G*^°^ is the initial Green's function used in the calculation of the collective 
modes, and {n-^ni) is the value obtained from spin and charge susceptibilities, then ^Tr [E'^^G*^"^] = U [n-^ri]) is 
satisfied exactly. The extent to which ^Tr [E'^^G'^^] (computed with G*^^^ instead of G*-*^^) differs from U {n-^ni) can 
then be used both as an internal accuracy check and as a way to improve the vertex corrections. 

If one is interested only in two-particle properties, namely spin and charge fluctuations, then this approach has 
the simple physical appeal of RPA but it satisfies key constraints that are always violated by RPA, namely the 
Mermin- Wagner theorem and the Pauli principle. To contrast it with usual RPA, that has a self-consistency only at 
the single-particle level, we call it the Two-Particle Self-Consistent approach (TPSC). [|| ^ [|| The TPSC gives a 
quantitative description of the Hubbard model not only far from phase transitions, but also upon entering the critical 
regime. Indeed we have shown quantitative agreement with Monte Carlo simulations of the nearest-neighbor p9| ] and 
next-nearest neighbor [p4| Hubbard model in two dimensions. Quantitative agreement is also obtained as one enters 
the narrow critical regime accessible in Monte Carlo simulations. We also have shown pS] in full generality that the 



TPSC approach gives the n ^ oo Umit of the O (n) model, while n = 3 is the physically correct (Heisenberg) limit. In 
two dimensions, we then recover both quantum-critical |]l9|] and renormalized classical |1§|] regimes to leading order in 
1/n. Since there is no arbitrariness in cutoff, given a microscopic Hubbard model no parameter is left undetermined. 
This allows us to go with the same theory from the non-critical to the beginning of the critical regime, thus providing 
quantitative estimates for the magnetic phase diagram of the Hubbard model, not only in two dimensions but also in 
higher dimensions |33|] . 

The main limitation of the approach presented in this paper is that it is valid only from weak to intermediate 
coupling. The strong-coupling case cannot be treated with frequency-independent irreducible vertices, as will become 
clear later. However, a suitable ansatz for these irreducible vertices in a Luttinger-Ward functional might allow us to 
apply our general scheme to this limit as well. 

Our approach predicts pO| that in two dimensions, Fermi liquid quasiparticles disappear in the renormalized classical 
regime £,afm oc exp{constjT), which always precedes the zero-temperature phase transition in two-dimensions. In this 
regime the antiferromagnetic correlation length becomes larger than the single-particle thermal de Broglie wave length 
£,th{— vf/T), leading to the destruction of Fermi Hquid quasiparticles with a concomitant appearance of precursors 
of antiferromagnetic bands ("shadow bands") with no quasi-particle peak between them. We stress the crucial role 
of the classical thermal spin fluctuations and low dimensionality for the existence of this effect and contrast our 
results with the earlier results of Kampf and Schrieffer ||35[| who used a susceptibility separable in momentum and 
frequency Xsp = /(q)5(^)- The latter form of Xsp = f (q) g(^) leads to an artifact that dispersive precursors of 
antiferromagnetic bands can exist at T = (for details see ^§1). We also contrast our results with those obtained in 
the fluctuation exchange approximation (FLEX), which includes self-consistency in the single particle propagators but 



neglects the corresponding vertex corrections. The latter approach predicts only the so-called "shadow feature" [g6 37 
which is an enhancement in the incoherent background of the spectral function due to antiferromagnetic fluctuations. 
However, it does not predict [Q the existence of "shadow bands" in the renormalized classical regime. These bands 
occur when the condition tu — ei^ ~ Ecr(k,w) -I- ^ = is satisfied. FLEX also predicts no pseudogap in the spectral 
function A{\iF,uj) at half-filling |Q. By analyzing temperature and size dependence of the Monte Carlo data and 
comparing them with the theoretical calculations, we argue that the Monte Carlo data supports our conclusion that 
the precursors of antiferromagnetic bands and the pseudogap do appear in the renormalized classical regime. We 
believe that the reason for which the FLEX approximation fails to reproduce this effect is essentially the same reason 
for which it fails to reproduce Hubbard bands in the strong coupling limit. More specifically, the failure is due to 
an inconsistent treatment of vertex corrections in the self-energy ansatz. Contrary to the electron-phonon case, these 
vertex corrections have a strong tendency to cancel the effects of using dressed propagators in the expression for the 
self-energy. 

Recently, there have been very exciting developments in photoemission studies of the High- Tc materials [ p9pO| ] that 
show the opening of the pseudogap in single particle spectra above the superconducting phase transition. At present, 
there is an intense debate about the physical origin of this phenomena and, in particular, whether it is of magnetic 
or of pairing origin. From the theoretical point of view there are a lot of formal similarities in the description of 
antiferromagnetism in repulsive models and superconductivity in attractive models. In Sec.(|V|) we use this formal 
analogy to obtain a simple analytical expressions for the self-energy in the regime dominated by critical pairing 
fiuctuations. We then point out on the similarities and differences in the spectral function in the case of magnetic 
and pairing pseudogaps. 

Our approach has been described in simple physical terms in Refs. p9] and [pOl. The plan of the present paper is 



as follows. After recalling the model and the notation, we present our theory in Sec. (III). There we point out which 
exact requirements of many-body theory are satisfied, and which are violated. Before Sec. (Ill), the reader is urged to 
read Appendix (Kh that contains a summary of sum rules, conservation laws and other exact constraints. Although 
this discussion contains many original results, it is not in the main text since the more expert reader can refer to the 
Appendix as need be. We also illustrate in this Appendix how an inconsistent treatment of the self-energy and vertex 
corrections can lead to the violation of a number of sum rules and inhibit the appearance of the Hubbard bands, a 



subject also treated in Sec.( Vl). Section ( |IVD compares the results of our approach and of other approaches to Monte 
Carlo simulations. We study in more details in Sec.(M) the renormalized classical regime at half-filling where, in 
two dimensions, Fermi liquid quasiparticles are destroyed and replaced by precursors of antiferromagnetic bands well 
before the T = phase transition. We also consider in this section the analogous phenomenon of pairing pseudogap 



which can appear in two dimensions when the pairing fluctuations become critical. The following section Sec. (VI) 



explains other attempts to obtain precursor s of antiferromagnetic bands and points out why approaches such as FLEX 



fail to see the effect. We conclude in Sec.(VH) with a discussion of the domain of validity of our approach and in 
Sec. (VIII) with a critical comparison with FLEX and pseudo-potential parquet approaches, listing the weaknesses and 
strengths of our approach compared with these. A more systematic description and critique of various many-body 
approaches, as well as proofs of some of our results, appear in appendices. 



II. MODEL AND DEFINITIONS 

We first present the model and various definitions. The Hubbard model is given by the Haniiltonian, 



H = - 2^ ijj y4aC]a + c]„c,„j +U2_^n^^n^l . (1) 

<ij><7 i 

In this expression, the operator Cia destroys an electron of spin a at site i. Its adjoint cj^ creates an electron and the 
number operator is defined by Uia = cli^Cia- The symmetric hopping matrix tij determines the band structure, which 
here can be arbitrary. Double occupation of a site costs an energy U due to the screened Coulomb interaction. We 
work in units where fc^ = 1, ?i = 1 and the lattice spacing is also unity, a — \. As an example that occurs later, the 
dispersion relation in the d-dimensional nearest-neighbor model is given by 

d 
ek = -2i^(cosfcj). (2) 

i=l 



A. Single-particle propagators, spectral weight and self-energy. 

We will use a "four"-vector notation k = (k,ifc„) for momentum-frequency space, and 1 = (ri,Ti) for position- 
imaginary time. For example, the definition of the single-particle Green's function can be written as 

G. (1,2) ^ - (r.ci. (n) 4 {t2)) ^ - {TrC (1) 4(2)) (3) 

where the brackets () represent a thermal average in the grand canonical ensemble, Tt is the time-ordering operator, 
and T is imaginary time. In zero external field and in the absence of the symmetry breaking G^ (1,2) — G^ (1—2) and 
the Fourier-Matsubara transforms of the Green's function are 

a(fc)=Ve-*--^ / dTe''^"^iG,(ri,ri)= /d(l)e-^'^(i)G,(l) (4) 

^-(i) = §E^^'^'^G'^w- (5) 

k 

As usual, experimentally observable retarded quantities are obtained from the Matsubara ones by analytical con- 
tinuation ikn —^ to -\- irj. In particular, the single-particle spectral weight A (k,a;) is related to the single-particle 
propagator by 

G.(k,.M=/|^^^ (6) 

J 2,11 ikn — LO 

A, (k,c^) = -2ImG? (k,c^) . (7) 

The self-energy obeys Dyson's equation, leading to 

ikn- (ek - Ai) - Sct (k,2fc„) 

It is convenient to use the following notation for real and imaginary parts of the analytically continued retarded 
self-energy 

E^ {k,ik„ ^Lo + n^) = ^'^ (k^Lu) + i^l (k,cj) . (9) 

Causality and positivity of the spectral weight imply that 



K iK^) < 0. 



(10) 



Finally, let us point out that for nearest- neighbor hopping, the Hamiltonian is particle-hole symmetric at hall-filling, 
Ckcr -^ cl+Qcr ! 4ia ^ Ck+Qff j wlth Q = (tt, tt) , implying that /j, = U/2 and that, 



G,(k,T)=-G,(k + Q,-r) 



(11) 



E(k,ifc„) - — 



S (k + Q, -iA:„) - ^ 



(12) 



B. Spin and charge correlation functions 



We shall be primarily concerned with spin and charge fluctuations, which are the most important collective modes 
in the repulsive Hubbard model. Let the charge and z components ol the spin operators at site i be given respectively 

by 



Pi (t) = n,^ (t) + n,i (r) 



(13) 



S^ = n^t (r) - n^l (t) , 



(14) 



The time evolution here is again that of the Heisenberg representation in imaginary time. 

The charge and spin susceptibilities in imaginary time are the responses to perturbations applied in imaginary-time. 
For example, the linear response of the spin to an external field that couples linearly to the z component 



is given by 

Xsp (y^i ^ Tj, Ti ^ Tj) = 

In an analogous manner, for charge we have 



^-^^^^^ = {TrSt in) S] (r,)) 



Xch \^i '^ji'^i ^j J 



HPj{tj)) 

min) 



{TTPt{Ti)pj{Tj)) -V? 



(15) 



(16) 



(17) 



Here n = {pi) is the filling so that the disconnected piece is denoted n^ . It is well known that when analytically 
continued, these susceptibilities give physical retarded and advanced response functions. In fact, the above two 
expressions are the imaginary-time version of the fluctuation-dissipation theorem. 

The expansion of the above functions in Matsubara frequencies uses even frequencies. Defining the subscript ch^ sp 
to mean either charge or spin, we have 



Xch,sp{(l,iqn) = 



duj' Xch,sp (q^^O 

TT Lu' - iq„ 



(18) 



X'4 (q,t) = i ([Pq it) , P-q (0)]) ; X% (q,0 = ^ ( [^q (t) ' "^-q (0)] ) ' 

The fact that Xch sp (q^^') i'' ^"^^l ^^^ odd in frequency in turn means that Xch,sp {<i,i<]n) is real 



Xch,spiq,iqn) = 



dw' ^'x'ch.sp (q^^O 

^ [io'f + iqnf 



(19) 



(20) 



a convenient feature for numerical calculations. The high-frequency expansion has 1/q^ as a leading term so that 
there is no discontinuity in Xch,sp (q, T)a.s. r ^ 0, contrary to the single-particle case. 



III. FORMAL DERIVATION 

To understand how to satisfy as well as possible the requirements imposed on many-body theory by exact results, 
such as those in Appendix (|A|), it is necessary to start from a general non-perturbative formulation of the many-body 
problem. We thus first present a general approach to many-body theory that is set in the framework introduced 
by Martin and Schwinger [Q, Luttinger and Ward |2^ and Kadanoff and Baym |g^ |2^]. This allows one to see 
clearly the structure of the general theory expressed in terms of the one-particle irreducible self-energy and of the 
particle-hole irreducible vertices. These quantities represent projected propagators and there is a great advantage in 
doing approximations for these quantities rather than directly on propagators. 

Our own approximation to the Hubbard model is then described in the subsection that follows the formalism. In 
our approach, the irreducible quantities are determined from various consistency requirements. The reader who is 
interested primarily in the results rather than in formal aspects of the theory can skip the next subsection and refer 
back later as needed. 

A. General formalism 

Following Kadanoff and Baym, ||2^ wc introduce the generating function for the Green's function 

In Z [0] = In ^r,e-^^(i)^^(2)^(T.2) ^ (21) 

where, as above, a bar over a number means summation over position and imaginary time and, similarly, a bar over a 
spin index means a sum over that spin index. The quantity Z is a functional of (f>a- , the position and imaginary-time 
dependent field. Z reduces to the usual partition function when the field ^o- vanishes. The one-particle Green's 
function in the presence of this external field is given by 

G.(1,2;[0]).-|MI^ (22) 

d<p„ (2,1) 

and, as shown by Kadanoff and Baym, the inverse Green's function is related to the self-energy through 

G^i = Gt7^-(/)-S. (23) 

The self-energy in this expression is a functional of 4>. 

Performing a Legendre transform on the generating functional In Z [0] in Eq. (UW with the help of the last two 
equations, one can find a functional <f> [G] of G that acts as a generating function for the self-energy 

The quantity $ [G] is the Luttinger- Ward functional, [p^ Formally, it is expressed as the sum of all connected skeleton 
diagrams, with appropriate counting factors. Conserving approximations start from a subset of all possible connected 
diagrams for $ [G] to generate both the self-energy and the irreducible vertices entering the integral equation obeyed 
by response functions. These response functions are then guaranteed to satisfy the conservation laws. They obey 
integral equations containing as irreducible vertices 

A complete and exact picture of one- and two-particle properties is obtained then as follows. First, the generalized 
susceptibilities Xa-a' (1, 3; 2) = —SG^ (1, 3) /6(f>a-' (2"*", 2) are calculated by taking the functional derivative of GG~^ and 
using the Dyson equation ( |23| ) to compute 5G~^ /Scf). One obtains |2^ 

X..'(1,3;2) = -G„(1,2)<5,^,,G,(2,3) + G, (l,2)r;V(2,3;4,5)x^,'(4,5;2)G, (3,3) (26) 

where one recognizes the Bethe-Salpeter equation for the three-point susceptibility in the particle-hole channel. The 
second equation that we need is automatically satisfied in an exact theory. It relates the self-energy to the response 
function just discussed through the equation 
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FIG. 1. The first line is a diagramatic representation of tlie Betlie-Salpeter equation Eq.(Pq) for tfie tfiree point susceptibility 
and the second line is the corresponding equation Eq.(E7|) for the self-energy. In the Hubbard model, the Fock contribution is 
absent, but in general it should be there. Solid lines are Green's functions and dashed lines represent the contact interaction 
U. The triangle is the three point vertex, while the three-point susceptibility x(l, 3; 2) is the triangle along with the attached 
Green's function. The usual two-point susceptibility is obtained by identifying points 1 and 3 in the Bethe-Salpeter equation. 
The rectangular box is the irreducible four-point vertex in the selected particle-hole channel. 

E,(l,2) = [/n_,5(l-2)-|-C/G,(l,2)r;V(2,2;4,5)x.'-.(4,5;l) (27) 

which is proven in Appendix (H). 

The diagrammatic representation of these two equations Eqs.(|2q)(p7|) appearing in Fig.(|l|) may make them look 
more familiar. Despite this diagrammatic representation, we stress that this is only for illustrative purposes. The rest 
of our discussion will not be diagrammatic. 

Because of the spin-rotational symmetry the above equations Eqs.(Efl) and (E^) can be decoupled into symmetric 
(charge) and antisymmetric (spin) parts, by introducing spin and charge irreducible vertices and generalized suscep- 
tibihties: 

-p — -p'ir j^ -pir . -p — -pir -r^ir /Oii\ 

1 c/i = -L II ~r i -|-| , 1 sp ^ i II 1 II \'^") 

Xch = '2{xTi+X„) ; Xsp = 2(x„ -XtJ (29) 

The usual two-point susceptibilities are obtained from the generalized ones as Xsp,c/i(l,2) = Xsp^cft(l, 1^; 2). The 
equation Eq.(|2^) for the generalized spin susceptibility leads to 

Xsp(l,3;2) = -2G(l,2)G(2,3)-r,p(2,3;4,5)G(l,2)G(3,3)x.p(4,5;2) (30) 

and similarly for charge, but with the plus sign in front of the second term. 

Finally, one can write the exact equation Eq. ( P7|) for the self-energy in terms of the response functions as 

E, (1,2) = Un^Jil -2) + J [r,p(2,2;4,5)x.p(4,5; 1) -H rch(2, 2; 4, 5)xc/.(4, 5; 1)] G, (l,2) . (31) 

Our two key equations are thus those for the three-point susceptibilities, Eq.(30), and for the self-energy, Eq.(Bl]). 
It is clear from the derivation in Appendix (|^) that these equations are intimately related. 

B. Approximations through local irreducible vertices. 



1. Conserving approximation for the collective modes. 

In formulating approximation methods for the many-body problem, it is preferable to confine our ignorance to 
high-order correlation functions whose detailed momentum and frequency dependence is not singular and whose 
influence on the low energy Physics comes only through averages over momentum and frequency. We do this here 
by parameterizing the Luttinger-Ward functional by two constants ri| and TV'^ . They play the role of particle- hole 
irreducible vertices that are eventually determined by enforcing sum rules and a self-consistency requirement at the 
two-particle level. In the present context, this functional can be also considered as the interacting part of a Landau 
functional. The ansatz is, 

•^ [G] = \g^ (T, T^) r^G^ (T, T^) +\g^ (T, T^) r?_^G_^ (T, T+) . (32) 

As in every conserving approximation, the self-energy and irreducible vertices are obtained from functional derivatives 
as in Eq.(|4|) and Eq.(p5|) and then the collective modes are computed from the Bethe-Salpeter equation Eq. (|30|) . The 
above Luttinger-Ward functional gives a momentum and frequency independent self-energy E^ , that can be absorbed 
in a chemical potential shift. From the Luttinger-Ward functional, one also obtains two local particle-hole irreducible 
vertices FJT^ and FJT^ 



crcr 



r-^(2,3;4,5)^^^=<S(2-5)<5(3-4)<5(4+-5)F-,. (33) 

We denote the corresponding local spin and charge irreducible vertices as 

[/.p EE F-_^ - F- ; [/,,^F-_,+F-^. (34) 

Notice now that there are only two equal-time, equal-point (i.e. local) two-particle correlation functions in this 
problem, namely {n-^rii) and (n?) = (n?) = {n^) = n/2. The last one is completely determined by the Pauli principle 
and by the known filling factor, while C/(n|nj) is the expectation value of the interaction term in the Hamiltonian. 
Only one of these two correlators, namely U{n^n \)^ is unknown. Assume for the moment that it is known. Then, we 



can use the two sum rules Eqs.(Al8) and (A17) that follow from the fluctuation-dissipation theorem and from the 
Pauli principle to determine the two trial irreducible vertices from the known value of this one key local correlation 
functions. In the present notation, these two sum rules are 

Xch (l, 1^) = ]^ X! X! XchiqMn) = {n-\) + {ni) + 2 (n|nx) - n^ (35) 

q iq,z 

Xsp (1, 1^) = ]^ XI XI Xspi^,iqn) = ("t) + ("i> - 2 (n^ni) (36) 

q iQr, 

and since the spin and charge susceptibilities entering these equations are obtained by solving the Bethe-Salpeter 
Eq.(pO[) with the constant irreducible vertices Eqs.(|33|)(p4|) we have one equation for each of the irreducible vertices 

n + 2{n,nCj -n^-^T. ^ J'^'^ ,y (37) 

We used our usual short-hand notation for wave vector and Matsubara frequency q — (q,i(7„). Since the self-energy cor- 
responding to our trial Luttinger-Ward functional is constant, the irreducible susceptibilities take their non- interacting 
value xo(g)- 



The local Pauli principle (n?) = (nj) leads to the following important sum-rule 



^ X X [^«P (^' *9") + 'X'^^ (^' **')] ^2n-r?, (39) 



q ig. 



which can be obtained by adding Eqs.(|38|),(37). This sum-rule imphes that effective interactions for spin Usp and 
charge Uch channels must be differe nt f rom one another and hence that ordinary RPA is inconsistent with the Pauli 
principle ffor details see Appendix ( |A 3| ) ) . 



Eqs.(p7|) and (38) determine Usp and Uch as a function of double occupancy (?i|nj^). Since double occupancy is a 
local quantity it depends little on the size of the system. It could be obtained reliably from a number of approaches, 
such as for example Monte Carlo simulations. However, there is a way to obtain double-occupancy self-consistently 
p9[ without input from outside of the present theory. It suffices to add to the above set of equations the relation 

Usp^gniO)U ; g,,(0)^^|^. (40) 

Eqs.(p8|) and (Em then define a set of self-consistent equations for U^p that involve only two-particle quantities. This 



ansatz is motivated by a similar approximation suggested by Singwi et al. |31| in the electron gas, which proved to be 
quite successful in that case. On a lattice we will use it for n < 1. The case n > 1 can be mapped on the latter case 
using particle-hole transformation. In the context of the Hubbard model with on-site repulsion, the physical meaning 
of Eq. (Eg) is that the effective interaction in the most singular spin channel, is reduced by the probability of having 
two electrons with the opposite spins on the same site. Consequently, the ansatz reproduces the Kanamori-Brueckner 
screening that inhibits ferromagnetism in the weak to intermediate coupling regime (see also below). We want to 
stress, however, that this ansatz is not a rigorous result like sum rules described above. The plausible derivation of 
this ansatz can be found in Refs. pl| , [ p9[ as well as, in the present notation, in Appendix (0). 

We have called this approach Two-Particle Self-Consistent to contrast it with other conserving approximations 
like Hartree-Fock or Fluctuation Exchange Approximation (FLEX) p3 that are self-consistent at the one-particle 
level, but not at the two-particle level. This approach |^ to the calculation of spin and charge fluctuations satisfles 
the Pauli principle (n^) = (rio-) = n/2 by construction, and it also satisfies the Mermin- Wagner theorem in two 
dimensions. 

To demonstrate that this theorem is satisfied, it suffices to show that {n^ni) — .gti (0) (^*T)("'i) does not grow 



indefinitely. (This guarantees that the constant C appearing in Eq.(A24) is finite.) To see how this occurs, write the 
self-consistency condition Eq.(|3l) in the form 



ritni) — — > ; ^ . 41) 



Consider increasing {n-^ni) on the right-hand side of this equation. This leads to a decrease of the same quantity on 
the left-hand side. There is thus negative feedback in this equation that will make the self-consistent solution finite. A 
more direct proof by contradiction has been given in Rcf. p9| : suppose that there is a phase transition, in other words 
suppose that (n|) (rij) = ^f/(fttftx)Xo(9)- Then the zero-Matsubara frequency contribution to the right-hand side of 
Eq.(Bl]) becomes infinite and positive in two dimensions as one can see from phase-space arguments (See Eq.(A24)). 



This implies that {n-^ni) on the left-hand side must become negative and infinite, but that contradicts the starting 
hypothesis since (n-f) (n|) = \U {nT^n\^)xo{q) means that (n-^ui) is positive. 

AlthoTjgh there is no finite-temperature phase transition, our theory shows that sufficiently close to half-filling (see 



Sec.( [VC )) there is a crossover temperature Tx below which the system enters the so-called renormali zed cl assical 



regime, where antiferromagnetic correlations grow exponentially. This will be discussed in detail in Sec.(VAl). 

Kanamori-Brueckner screening is also included as we already mentioned above. To see how the screening occurs, 
consider a case away from half-filling, where one is far from a phase transition. In this case, the denominator in the 
self-consistency condition can be expanded to linear order in U and one obtains 

("t) w) l + AC/ 



where 



^-i^T^Xoicir- (43) 



n^ N 

q 



Clearly, quantum fluctuations contribute to the sum appearing above and hence to the renormalization of Usp — 
9U (0) ^- ^^^ value of A is found to be near 0.2 as in explicit numerical calculations of the maximally crossed 
Kanamori-Brueckner diagrams. [Q At large U, the value of Usp = g-\i (0) U ^ 1/A saturates to a value of the order 
of the inverse bandwidth which corresponds to the energy cost for creating a node in the two-body wave function, in 
agreement with the Physics described by Kanamori. ||] 
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FIG. 2. Dependence on U of the charge and spin effective interactions (irreducible vertices). The temperature is chosen 
so that for aU U , it is above the crossover temperature. In this case, temperature dependence is not significant. The filling is 
n= 1 

To illustrate the dependence of C/jp, Uck on bare U we give in Fig.(0) a plot of these quantities at half- filling where 
the correlation effects are strongest. The temperature for this plot is chosen to be above the crossover temperature Tx 
to the renormalized classical regime, in which case the dependence of \J sp and Vch on temperature is not significant. 
As one can see, t/^p rapidly saturates to a fraction of the bandwidth, while Vch rapidly increases with f/, reflecting the 
tendency to the Mott transition. We have also shown previously in Fig. (2) of Ref. [g9| that Vg-p depends only weakly 
on filling. Since C/^p saturates as a function of \J due to Kanamori-Brueckner screening, the crossover temperature 
Tx also saturates as a function of U . This is illustrated in Fig.(Q) along with the mean- field transition temperature 
that, by contrast, increases rapidly with U . 

Quantitative agreement with Monte Carlo simulations on the nearest-neighbor p9| and next-nearest-neighbor mod- 
els p3] is obtained p^] for all fillings and temperatures in the weak to intermediate coupling regime U < 8t. This is 
discussed further below in Sec.(prv|). We have also shown that the above approach reproduces both quantum-critical 
and renormalized-classical regimes in two dimensions to leading order in the 1/n expansion (spherical model) p3]. 

As judged by comparisons with Monte Carlo simulations |4q|, the particle-particle channel in the repulsive two- 
dimensional Hubbard model is relatively well described by more standard perturbative approaches. Although our 
approach can be extended to this channel as well, we do not consider it directly in this paper. It manifests itself only 
indirectly through the renormalization of Ugp and Uch that it produces. 



2. Single-particle properties 



As in any implementation of conserving approximations, the initial guess for the self-energy, E^'^\ obtained from 
the trial Luttinger-Ward functional is inconsistent with the exact self-energy formula Eg. (|3l|) . The latter formula 
takes into account the feedback of the spin and charge collective modes actually calculated from the conserving 
approximation. In our approach, we use this self-energy formula Eq.(|3l|) in an iterative manner to improve on 
our initial guess of the self-energy. The resulting formula for an improved self-energy E'^' has the simple physical 
interpretation of paramagnon theories. |46 

As another way of Physically explaining this point of view, consider the following: The bosonic collective modes are 
weakly dependent on the precise form of the single-particle excitations, as long as they have a quasiparticle structure. 
In other words, zero-sound or paramagnons exist, whether the Bethe-Salpeter equation is solved with non-interacting 
particles or with quasiparticles. The details of the single-particle self-energy by contrast can be strongly influenced 
by scattering from collective modes because these bosonic modes are low-lying excitations. Hence, we first compute 
the two-particle propagators with Hartree-Fock single-particle Green's functions, and then we improve on the self- 
energy by including the effect of collective modes on single-particle properties. The fact that collective modes can be 
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FIG. 3. Crossover temperature at half-filling as function of U compared with the mean-field transition temperature 



calculated first and self-energy afterwards is reminiscent of renormalization group approaches, |47| B where collective 
modes are obtained at one-loop order while the non-trivial self-energy comes out only at two-loop order. 

The derivation of the general self-energy formula Eq. (|3^) given in Appendix (^) shows that it basically comes from 
the definition of the self-energy and from the equation for the collective modes Eq.(|30|). This also stands out clearly 
from the diagrammatic representation in Fig.(||). By construction, these two equations Eqs.(30) and ( pl| ) satisfy the 
consistency requirement -^TrYiG = U (n^ni) (see Appendix (|^)), which in momentum and frequency space can be 
written as 



k 



{k)G^{k)e-"'"^ = U{n^ni). 



(44) 



The importance of the latter sum rule, or consistency requirement, for approximate theories should be clear from 
the appearance of the correlation function (nini) that played such an important role in determining the irreducible 
vertices and in obtaining the collective modes. Using the fiuctuation dissipation theorem Eqs.(|3q),(|35|) this sum-rule 
can be written in form that explicitly shows the relation between the self-energy and the spin and charge susceptibilities 



^ J2 i^Ak) - Un_.] G.(fc) = 7^ E [^-'^(«) - ^«p(9)] 



N 



(45) 



To keep as much as possible of this consistency, we use on the right-hand side of the self-energy expression Eq. (|3l| ) 
the same irreducible vertices and Green's functions as those that appear in the collective- mode calculation Eq. (pOP . 
Let us call G*^"^ the initial Green's function corresponding to the initial Luttinger-Ward self-energy E^°). Our new 
approximation for the self-energy E*-^' then takes the form 



U T 



41) (fc) - C/n_. + Tl^ E PspXspiq) + UchXchiq)] G^^\k + q). 



(46) 



Note that S^ {k) satisfies particle- hole symmetry Eq.([l^) where appropriate. This self-energy expression ( p^ ) is 
physically appealing since, as expected from general skeleton diagrams, one of the vertices is the bare one [/, while 
the other vertex is dressed and given by Ugp or Uch depending on the type of fluctuation being exchanged. It is 
because Migdal's theorem does not apply for this problem that Ugp and Uch are different from the bare U at one of 
the vertices. Usp and Uch bere take care of vertex corrections. [ p8[ 

The use of the full Ga{k -\- q) instead of G"(fc + q) in the above expression Eq.(p6|) would be inconsistent with 
frequency- independent irreducible vertices. For the collective mode Eq.(30) this is well known to lead to the violation 
of the conservation laws as was discussed in detail in the previous subsection. Here we insist that the same is true in 
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the calculation of the effect of electronic collective modes on the single- particle properties. Formally, this is suggested 
by the similarity between the equation for the susceptibility Eq.(|3y) and that for the self-energy Eq.(|3l|) in terms of 
irreducible vertices. More importantly, two physical effects would be absent if one were to use full G and frequency 
independent irreducible vertices. First, upper and lower Hubbard bands would not appear because the U'^ /lu high- 
frequency behavior in Eg . (pSJ) that is necessary to obtain the Hubbard bands would set in too late, as we discuss in 



Sec.(A2) and in Sec. (VIA). This result is also apparent from the fact that FLEX calculations in infinite dimension 
do not find upper and lower Hubbard bands pt9| ] where the exact numerical solution does. The other physical effect 
that would be absent is precursors of antiferromagnetic bands, Sec.(V| and the pseudogap in yl(kF,w), that would 



not appear for reasons discussed in Sec.(V]). We also will see in Sec. (IV) below that FLEX calculations of the single- 
particle Green's function, significantly disagree with Monte Carlo data, even away from half-filling, as was already 
shown in Fig. 1 of Ref. |§. 

The chemical potential for interacting electrons fi is found from the usual condition on particle number 

n = ^j:Gi^Hk)eM~^k^O-) -^E-- ^"P(-f"»') (47) 

This chemical potential // is, of course, different from fig but the Luttingcr sum rule X]^(~^k + /^ ^ S^^-*) = ria 
is satisfied to a high accuracy (about few percent) for all fillings and temperatures Tx < T <C W. As usual this 
occurs because the change in /i*^^^ — hq is compensated by the self-energy shift on the Fermi surface Ti^^^kp, 0). For 
T < Tx there is some deviation from the Luttinger sum rule which is due to the appearance of the precursors of the 
antiferromagnetic bands below Tx (Sec. 0) which develop into true SDW bands at T = 0. 

It is important to realize that G^"' on the right hand side of the equation for the self-energy S cannot be calculated 
as G^"^ — l/{uj — £k + Z^^^-'), because otherwise it would not reduce to zero-temperature perturbation theory when it 



is appropriate. As was pointed out by Luttinger, (see also section A 4) the "non-interacting" Green's function used in 
the calculation for E should be calculated as G*^°^ = l/(^ — Ek — S'-"-*(ki?,0) -I- A^^"-*), where /j,*^"^ is calculated on the 
same level of accuracy as T,^'^\]s.f,0), i.e. from Eq.(|47|) with E'^"^(k, i/c„). In our calculation below, we approximate 
^(1) _ E(^^(ki?,0) by fiQ because for the couphng strength and temperatures considered in this paper {U < W/2, 
Tx < T ^ W) the Luttinger theorem is satisfied to high accuracy and the change of the Fermi surface shape is 
insignificant. In addition, at half-filling the condition /i — S(kj?, 0) = /xg is satisfied exactly at any U and T because of 
particle-hole symmetry. For somewhat larger coupling strengths and away from half-filling, one may try to improve 
the theory by using G^"' = l/(cj — £k — 5](-^)(kF,0) -I- ^^^^), with E^^^ and fi found self-consistently. However, the 
domain of validity of our approach is limited to the weak-to-intermediate coupling regime since the strong-coupling 
regime requires frequency-dependent pseudopotentials (see below). 

Finally, let us note that, in the same spirit as Landau theory, the only vertices entering our theory are of the type 
Fij^ and Fn, or, through Eq.(p4), Ugp and Uch- In other words, we look at the problem from the longitudinal spin 
and charge particle-hole channel. Consequently, in the contact pseudopotential approximation the exact equation 
for the self-energy Eq.(pl|) reduces to our expression Eq. (|4^) which does not have the factor 3 in the front of the 



spin susceptibility. This is different from some paramagnon theories, in which such factor was introduced to take 



care of rotational invariance. However, we show in Appendix (El) that these paramagnon theories are inconsistent 
with the sum-rule Eq.([45|) which relates one and two-particle properties. In our approach, questions about transverse 
spin fluctuations are answered by invoking rotational invariance xfp = xip — xtp- I'^ particular, one can write the 



expression for the self-energy Eq.(46) in an explicitly rotationally invariant form by replacing Xsp by {l/S)Tr[x'^'^ 



sp J 



If calculations had been done in the transverse channel, it would have been crucial to do them while simultaneously 
enforcing the Pauli principle in that channel. In functional integration methods, it is well known that methods that 
enforce rotational invariance without enforcing the Pauli principle at the same time give unphysical answers, such as 
the wrong factor 2/3 in the RPA susceptibility pJ Xsp — X^/{^ ~ (2/3)t/x'') or wrong Hartree-Fock ground state. 



3. Internal accuracy check 

The quantitative accuracy of the theory will be discussed in detail when we compare with Monte Carlo calculations 
in the next section. Here we show that we can use the consistency requirement between one- and two-particle properties 
Eq. ( [44| ) to gauge the accuracy of the theory from within the theory itself. 

The important advantage of the expression for the self-energy So- (k) given by Eq.(|4q) is that, as shown in Ap- 
pendix (H) , it satisfies the consistency requirement between one- and two-particle properties Eq. (E3) , in the following 
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r_o- TV ^-^ 

k 



S^i) (k) Gi°) (fc) e-*'^"" = 



C/ 



>'^T'^i> 



(48) 



Let Ga be defined by 






= Gq ^ — E*^^). We can use the fact that in an exact theory we should have 



Tr 



y(l)/^(l) 



in the above expression instead of Tr 



compute by how much Tr 






differs from Tr 



v(i)r^(o) 



to check the accuracy of the theory. It suffices to 
. In the parameter range U < At and n,T arbitrary 



but not too deep in the, soon to be described, renormalized-classical regime, we find that Tr 



v(i)r^(o) 



differs from 



Tr 



v(i)r^(i) 



by at most 15%. Another way to check the accuracy of our approach is to evaluate the right-hand side 



.(1) 



of the /—sum rule Eqs.(A25) with nk^ — Ga (k, 0^) and to compare with the result that had been obtained with 
/k.CT- Again we find the same 15% disagreement, at worse, in the same parameter range. As one can expect, this 
deviation is maximal at half- filling and becomes smaller away from it. 

Eq.(p6|) for the self-energy E^^^ already gives good agreement with Monte Carlo data but the accuracy can be 

improved even further by using the general consistency condition Eq. (B3) on Tr Sjj Go- to improve on the ap- 
proximation for vertex corrections. To do so we replace Ugp and UcH on the right-hand side of Eq.(Effl by aUsp and 

aUch with a determined self-consistently in such a way that Eq.(|48|) is satisfied with Ga (k) replaced by Ga (k). 
For C/ < 4, we have a < 1.15. The slight difference between the irreducible vertices entering the collective modes 
and the vertex corrections entering the self-energy formula can be understood from the fact that the replacement of 
irreducible vertices by constants is, in a way, justified by the mean-value theorem for integrals. Since the averages 
are not taken over the same variables, it is clear that the vertex corrections in the self-energy formula and irreducible 
vertices in the collective modes do not need to be strictly identical when they are approximated by constants. 

Before we move on to comparisons with Monte Carlo simulations, we stress that E^^^ given by Eq. (Eo) cannot be 
substituted back into the calculation of Xsp,ch by simply replacing xo = GqGq with the dressed bubble xo — GG. 
Indeed, this would violate conservation of spin and ch arge and /-sum rule. In particular, the co ndit ion Xsp,c/i(q = 
0,zg„ f^ 0) = that follows from the Ward identity ( A31 ) would be violated as we see in Eq.( A26 ). In the next 
order, one is forced to work with frequency-dependent irreducible vertices that offset the unphysical behavior of xq 
at non-zero frequencies. 



IV. NUMERICAL RESULTS AND COMPARISONS WITH MONTE CARLO SIMULATIONS 

In this section, we present a few numerical results and comparisons with Monte Carlo simulations. We divide this 
section in two parts. In the first one we discuss data sufficiently far from half-filling, or at high enough temperature, 
where size effects are unimportant for systems sizes available in Monte Carlo simulations. In the second part, we discuss 
data at half-filling. There, size effects become important below the crossover temperature Tx where correlations start 
to grow exponentially (Sec.(0)). All single-particle properties are calculated with our approximation Eq.(^ for the 
self-energy using the vertex renormalization a explained in the previous section. The results would differ at worse by 
15% if we had used a = 1. 



A. Far from the crossover temperature Tx 



1. Two-particle properties 



We have shown previously in Fig.4a-d of Ref. p3] and in Figs. 2-4 and Fig. 6 of Ref. pi] that both spin and charge 
structure factor sufficiently away from the crossover temperature Tx are in quantitative agreement with Monte Carlo 
data for values of U as large as the bandwidth. On the qualitative level, the decrease in charge fluctuations as one 
approaches half-filling has been explained p9[ | as a consequence of the Pauli principle embodied in the calculation of 
the irreducible vertex Uch- ^^ 

Here we present on Fig.(0) and Fig.(p|) comparisons with a dynamical quantity, namely the spin susceptibility. 
Similar comparisons, but with a phenomenological value of C/sp, have been done by Bulut et al. Ref. |52|. Fig.(U) 
shows the staggered spin susceptibility as a function of Matsubara frequencies for n = 0.87, T = 0.25 and U — A. The 
effect of interactions is already quite large for the zero-frequency susceptibility. It is enhanced by a factor of over 5 
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FIG. 4. Comparisons between Monte Carlo simulations [M and our theory for the spin susceptibility at Q = (tt, tt) as a 
function of Matsubara frequency. The temperature is T = 0.25, and the system size 8x8. The factor 1/2 on the vertical axis 
is due to the fact that the susceptibility in MM is xh a quantity that is by definition twice smaller then ours and that of |53] 



compared with the non-interacting value. Nevertheless, one can see that the theory and Monte Carlo simulations are 
in good agreement. 

Fig.(m shows the temperature dependence of the zero- frequency staggered spin susceptibility for the same filling 
and interaction as in the previous figure. Symbols represent Monte Carlo simulations from Refs. [Q and |5^], the solid 
line is for our theory while dotted and dashed lines are for two versions of FLEX. Surprisingly, the fully conserving 
FLEX theory, (dashed line) compares worse with Monte Carlo data than the non-conserving version of this theory 
that neglects the so-called Aslamasov-Larkin diagrams (dotted line). By contrast, our theory is in better agreement 
with the Monte Carlo data than FLEX for the staggered susceptibility Xsp (q = ("""i ti") ,i^n — 0), and at the same 
time it agrees exactly with the conservation law that states that Xsp,ch (q = 0,iw,i ^ 0) = 0. 

Finally, Fig.(||) shows the double occupancy {n-^rii) as a function of filling for various values of U. The symbols 
again represent Monte Carlo data for T = 1/6, and the lines are the results of our theory. Everywhere the agreement 
is very good, except for n = 1,C/ = 4. In the latter case, t he sy stem is already below the crossover temperature 
Tx to the renormalized classical regime. As explained in Sec. (VI]), the appropriate procedure for calculating double 
occupancy in this case is to take for (n^ni) its value (dotted line) at Tx instead of using the ansatz Eq.(E^). In any 
case, the difference is not large. 



2. Single-particle properties 



Fig. 1(a) of Ref. pfl| shows G(k,r) for filling n — 0.875, temperature T = 0.25 and U — 4 for the wave vector on 
the 8x8 lattice which is closest to the Fermi surface, namely (tt, 0). Our theory is in agreement with Monte Carlo 
data and with the parquet approach p3[ but in this regime second-order perturbation theory for the self-energy gives 
the same result. Surprisingly, FLEX is the oiily theory that disagrees significantly with Monte Carlo data. The good 
performance of perturbation theory (see also |p4| ) can be explained in part by compensation between the renormalized 
vertices and susceptibilities {Usp < U, Xspiq) > Xo{q)] Uch > U, Xch{q) < Xoiq) )■ 

We have also calculated Re(S (ifc„) /i/c„) and compared with the Monte Carlo data in Fig. 2a of Ref. ||5^] obtained 
at n = 0.87, U = 4, (3 = 6. Our approach agrees with Monte Carlo data for all frequencies, but again second-order 
perturbation theory gives similar results. 



B. Close to crossover temperature Tx at half- filling 
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FIG. 5. Comparisons between the Monte Carlo simulations (BW) and FLEX calculations presented in Fig. 19 of Ref. |53] 
and our theory for the spin susceptibility at Q = (n, n) as a function of temperature at zero Matsubara frequency. The filled 
circles (BWS) are from Ref. M 




FIG. 6. Comparisons between the Monte Carlo simulations of Ref. p^ and our theory (solid lines) for the filling dependence 
of the double occupancy. The results are are for T = 1/6 as a function of filling and for various values of U except for [7 = 4 
where the dotted line shows the results of our theory at the crossover temperature T — Tx ■ 
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1. Two-particle properties 

The occurrence of the crossover temperature Tx at half-filhng is perhaps best illustrated in the upper part of Fig. (|7|) 
by the behavior of the static structure factor Ssp {tt, tt) for [/ = 4 as a function of temperature. When the correlation 
length becomes comparable to the size of the system used in Monte Carlo simulations, [Q the static structure factor 
starts to increase rapidly, saturating to a value that increases with system size. The solid line is calculated from our 
theory for an infinite lattice. The Monte Carlo data follow our theoretical curve (solid line) until they saturate to 
a size-dependent value. The theory correctly describes the static structure factor not only above Tx bu t also as we 



enter the renormalized classical regime at Tx- Analytical results for this regime are given in Sec.(V_A_l). Note that 
the RPA mean-field transition temperature for this value of U is more than three times larger than Tx ~ 0.2. The 
size-dependence of Monte Carlo data for Ssp (q) at all other values of q 7^ (tt, tt) available in simulations is negligible 
and our calculation for infinite system reproduces this data (not shown). 

2. Single-particle properties 

Equal-time (frequency integrated) single-particle properties are much less sensitive to precursor effects than dynam- 
ical quantities as we now proceed to show. For example, n (k) — G (k, 0~) is a sum of G (k, ikn) over all Matsubara 
frequencies. We have verified (figure not shown) that j^ X^ko- "-ko-c^^ek/i^fcx obtained from Monte Carlo simulations 
p6[ is given quite accurately by either second-order perturbation theory or by our theory. This has very important 
consequences since, for this quantity, the non-interacting value differs from second-order perturbation theory by at 
most 15%. This means that the numerical value of the right-hand side of the / sum-rule Eq.(A25|) is quite close to 



that obtained from the left-hand side using our expression for the spin and charge susceptibility. 

One can also look in more details at n (k) itself instead of focusing on a sum rule. Fig.(@) shows a comparison of 
our theory and of second order perturbation theory with Monte Carlo data for n (k) obtained for a set of lattice sizes 
from 6 X 6 to 16 X 16 at n = 1, r = 1/6, U — 4. Size effects appear unimportant for this quantity at this temperature. 
These Monte Carlo data have been used in the past [Q to extract a gap by comparison with mean field SDW theory. 
Our theory for the same set of lattice sizes is in excellent agreement with Monte Carlo data and predicts a pseudogap 
at this temperature, as we will discuss below. However, for available values of k on finite lattices, second order 
perturbation theory is also in reasonable agreement with Monte Carlo data for n (k) . Since second order perturbation 
theory does not predict a pseudogap, this means that n (k) is not really sensitive to the opening of a pseudogap. This 
is so both because of the finite temperature and because the wave vectors closest to the Fermi surface are actually 
quite far on the appropriate scale. For this filling, the value of n (k) is fixed to 1/2 on the Fermi surface itself. 

It is thus necessary to find a dynamical quantity defined on the Fermi surface whose temperature dependence will 
allow us to unambiguously identify the pseudogap regime in both theory and in Monte Carlo data. The most dramatic 
effect is illustrated in the lower part of Fig.(0) where we plot the quantity z (T) defined by [Q Q 

The physical meaning of this quantity z (T) is that it is an average of the single-particle spectral weight A (lip, ^) within 
T = l/f3 of the Fermi level {uj — 0). When quasiparticles exist, this is the best estimate of the usual zero-temperature 
quasiparticle renormalization factor z = 1/(1 — dTi/duj) that can be obtained directly from imaginary-time Monte 
Carlo data. For non-interacting particles z (T) is unity. For a normal Fermi liquid it becomes equal to a constant less 
than unity as the temperature decreases since the width of the quasiparticle peak scales as T^ and hence lies within T 
of the Fermi level. However, contrary to the usual z = 1/(1 — dYi/duj) this quantity gives an estimate of the spectral 
weight A(ki?,u^ around the Fermi level, even if quasiparticles disappear and a pseudogap forms, as in the present 
case, (see Sec.(0)). 

One can clearly see from the lower part of Fig.(n) that while second-order perturbation theory exhibits typical 
Fermi-liquid behavior for z (T) , both Monte Carlo data (53| and a numerical evaluation of our expression for the 
self-energy lead to a rapid fall-off of z (T) below Tx (for [/ = 4, Tx « 0.2 Q). The rapid decrease of z (T) clearly 
suggests non Fermi-liquid behavior. We checked also that our theory reproduces the Monte Carlo size-dependence. 
This dependence is explained analytically in Sec^ V A 2|) . In Ref. [B0[ we have shown that at half-filling, our theory 



gives better agreement with Monte Carlo data |53(| for G (kf , r) than FLEX, parquet or second order perturbation 
theory. 

To gain a qualitative insight into the meaning of this drop in z (T) , we use the analytical results of the next section to 
plot in Fig.(|^) the value of A (kp, co) . This plot is obtained by retaining only the contribution of classical fluctuations 
Eq.(p9[) to the self-energy. One sees that above Tx, there is a quasiparticle but that at T ^ Tx a minimum instead 
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FIG. 7. The upper part of the figure, adapted from Ref. [p9|, shows the temperature dependence of S^^tTj-k) at half-fiUing 
n — 1. The sohd line is our theory for an infinite system while symbols are Monte Carlo data from Ref. |56|. The bottom part 
of the figure, adapted from Re f. [pC| , shows the behavior of z{T) — — 2G(kF,/3/2) in Eq.([49), as a function of temperature 
as obtained from Monte Carlo psfsimulations (symbols), from second order pertrubation theory (dashed line) and from our 
theory for an infinite system (solid line) and for a 16 x 16 lattice (dashed line). 
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FIG. 8. Occupation number n{k) as a function of wave vector k at half-filling for T = 1/6, U — 4, and system sizes 6x6 
to 16 X 16. The symbols are Monte Carlo results from Ref. [pTl while the solid line is our theory and the dotted line is the 
prediction from second order perturbation theory. The dashed line shows the result for L*' = as a reference. 

of a maximum starts to develop at the Fermi surface lo — 0. Below Tx, the quasiparticle maximum is replaced by two 
peaks that are the precursors of antiferromagnetic bands. This is discussed in detail in much of the rest of this paper. 



C. Phase diagram 



The main features predicted by our approach for the magnetic phase diagram of the nearest-neighbor hopping 
model have been given in Ref. [p3 . Needless to say, all our considerations apply in the weak to intermediate coupling 
regime. Note also that both quantum critical and renormalized classical properties of this model have been studied 
in another publication 133]. The shape of the phase diagram that we find is illustrated in Fig.(nO) for U = 2.5 and 
[/ = 4. 

At zero temperature and small filling, the system is a paramagnetic Fermi liquid, whatever the value of the interac- 
tion U (< W). Then, as one moves closer to half-filling, one hits a quantum critical point at a value of filling ric- Since, 
Usp in our theory saturates with increasing U, the value of ric is necessarily larger than about nc{U = oo) = 0.68. 
At this point, incommensurate order sets in at a wave vector {qc,Tr) or at symmetry-related points. Whatever the 
value of U, the value of qc is contained [p^ in the interval 0.747r < Qc < tt, increasing monotonously towards O.Tin 
as U increases. Since our approach applies only in the paramagnetic phase, at zero temperature we cannot move 
closer to half-filling. Starting from finite-temperature then, the existence of long-range order at low temperature is 
signaled by the existence of a crossover temperature Tx {n, U) below which correlations start to grow exponentially. 
We have already discussed the meaning of Tx {n, U) at half-filling. This crossover temperature becomes smaller and 
smaller as one moves away from half-filling, until it reaches the quantum-critical point that we just discussed. The 
correlations that start to grow at Tx {n, U) when n = 1 are at the antiferromagnetic wave vector, and they stay at 
this wave vector for a range of fillings n. Finally, at some filling, the correlations that start to grow at Tx {n, U) are 
at an incommensurate value until the quantum-critical point is reached. 

Note that the above phase diagram is quite different from the predictions of Hartree-Fock theory mostly because 
of the strong renormalization of Usp- This quantitative change leads to qualitative changes in the Hartree-Fock phase 
diagram since, for example, Stoner ferromagnetism never occurs in our picture. While the existence of ferromagnetism 
in the strong coupling limit has been proven only recently [ p9| , the absence of Stoner ferromagnetism in the Hubbard 
model was already suggested by Kanamori a long time ago and was verified by more recent studies. |44| |6y] |6^] 
More relevant to the present debate though, is the fact that SDW order persists away from half-filling for a finite 
range of dopings. While this is in agreement with slave-boson approaches |62|] and studies |63| ] using the infinite- 
dimension methodology, |ll| it is in clear disagreement with Monte Carlo simulations |64|. Our approach certainly 
fails sufficiently below Tx, but given the successes described above, we believe that it can correctly predict the 
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FIG. 9. Qualitative sltetch of the spectral weight at the Fermi wave vector at half-filling for three temperatures. This plot 
is obtained by retainig only the classical contribution to the self-energy using parameters corresponding to the typical U — 4 
of Monte Carlo simulations. The top plot is for T > Tx, the middle one for T ~ Tx and the bottom one for T < Tx- The 
precursors of antiferromagnetic bands would look like this last figure. 
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FIG. 10. Crossover temperature Tx as a function of filling for t/ = 4 and U = 2.5. On this crossover line, ^ is enhanced by 
a factor of 500 over the bare value. Filled symbols indicate that the crossover is at the antiferromagnetic wa ve ve ctor, while 
open symbols indicate a crossover at an incommensurate wave vector. Reproduced with permission from Ref. (iM 
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exponential growth of fluctuations at Tx ■ It would be difficult to imagine how one could modify the theory in such 
a way that the growth of magnetic fluctuations does not occur even at incommensurate wave vectors. Also, such an 
approach would also need to stop the growth of fluctuations that we find as we approach the quantum critical point 
along the zero temperature axis, from the low-filling, paramagnetic side, where Tx (ji < ric, U) = 0. 

It could be that Monte Carlo simulations |64| fail to see long-range order at zero temperature away from half-filling 
because at zero temperature, in the nearest-neighbor model, this order has a tendency to being incommensurate 
everywhere except at n = 1. Furthermore, as we saw above, this incommensuration is in general far from one of 
the available wave vectors on an 8 x 8 lattice. It comes close to (O.TStt, tt) only for the largest values of U available 
by Monte Carlo. Hence, incommensurate order on small lattices is violently frustrated not only by the boundary 
conditions, but also by the fact that there is no wave vector on what would be the Fermi surface of the infinite system. 
This means that the electron-electron interaction scatters the electrons at wave-vectors that are not those where the 
instability would show up, rendering these scatterings not singular. This is clearly an open problem. 

V. REPLACEMENT OF FERMI LIQUID QUASIPARTICLES BY A PSEUDOGAP IN TWO 

DIMENSIONS BELOW Tx- 

One of the most striking consequences of the results discussed in the context of Monte Carlo simulations is the fall 
of the spectral weight below the temperature Tx where antiferromagnetic fluctuations start to grow exponentially in 
two dimensions. We have already shown in a previous publication [ pO[ that this corresponds to the disappearance 
of Fermi liquid quasiparticles at the Fermi surface, well above the zero temperature phase transition. We also found 
that, simultaneously, precursors of the antiferromagnetic bands develop in the single-particle spectrum. Given the 
simplicity of our approach, it is possible to demonstrate this phenomenon analytically. This is particularly important 
here because size effects and statistical errors make numerical continuation of the Monte Carlo data to real frequencies 
particularly difficult. Such analytic continuations using the maximum entropy method JM have, in the past, lead to 
a conclusion different from the one obtained later using singular value decomposition ||65|| 7 

In this section then, we will consider the conditions for which Fermi liquid quasiparticles can be destroyed and 
replaced by a pseudogap in two dimensions. The major part of this section will be concerned with the single particle 
pseudogap and the precursors of antiferromagnetic bands in the vicinity of the zero temperature antiferromagnetic 
phase transition in the positive U Hubbard model. However, it is well known that the problem of superconductivity is 
formally related to the problem of antiferromagnetism, in particular at half-filling where the nearest-neighbor hopping 
positive \J Hubbard model maps exactly onto the nearest-neighbor negative U Hubbard model. The corresponding 
canonical transformation maps the q = (tt, tt) spin correlations of the repulsive model onto the q = pairing and 
q = (tt, tt) charge correlations of the attractive model while the single-particle Green's functions of both models are 
identical. Thus all our results below concerning the opening of the pseudogap in A (k^?, w) in the repulsive U half- filled 
Hubbard model are directly applicable to the attractive V model at half-filling, the only difference being in the physical 
interpretation. While in the case of repulsive interaction the pseudogap is due to the critical thermal spin fluctuation, 
in the case of attractive interactions it is, obviously, due to the critical thermal pairing and charge fluctuations. Away 
from half-filling the mapping between two models is more complicated and the single particle spectra in the pairing 
pseudogap regime Ai^p^u:^ have important qualitative differences with the single particle spectra in the magnetic 
pseudogap regime. However, even in this case there are very useful formal similarities between two problems so that 



in Subsec.(VF) we will give some simple analytical results for the self-energy in the regime dominated by critical 
pairing fluctuations. 

The problem of precursor effects in the repulsive Hubbard model has been first studied by Kampf and Schrieffer 
p5[ . Their analysis however was done at zero temperature and although the precursor effect that they found, called 
"shadow bands" , looks similar to what we find, there are a number of important differences. For example, they 
find a quasiparticle between the precursors of antiferromagnetic bands, while we do not. Also, one does not obtain 
precursors at zero temperature when one uses our more standard expression for the dynamical susceptibility instead of 
the phenomenological form XKShr = /(ci)ff(^) that they use. The physical reason why a function that is separable in 
both momentum and frequency, such as XKShr, leads to qualitatively different results than the conventional one has 
been explained in Ref. ]3q |. The microscopic justification for XKShr is unclear. We comment below on this problem 
as well as on some of the large related literature that has appeared lately. 

Repeating some of the arguments of Ref. Q , we first show by general phase space arguments that the feedback of 
antiferromagnetic fluctuations on quasiparticles has the potential of being strong enough to destroy the Fermi liquid 
only in low enough dimension, the upper critical dimension being three. Then we go into more detailed analysis 
to give explicit analytic expressions for the quasi-singular part of the self-energy, first in Matsubara frequency. The 
analysis of the self-energy expression directly in real- frequencies is in Appendix (O) . The latter analysis is useful to 
exhibit in the same formalism both the Fermi liquid limit and the non- Fermi liquid limit. 
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For simplicity we give asymptotics for n = 1 at the Fermi wave vector, where eCkp) — 0, but similar results apply 
for n ^ 1 as long as there is long-range order at T = and one is below Tx- This case is also discussed briefly, but 
for more details the reader is referred to Ref. |Q . 

A. Upper critical dimension for the destruction of quasiparticles by critical fluctuations. 

Before describing the effect of spin fluctuations on quasiparticles, we first describe the so-called renormalized classical 
regime of spin fluctuations that precedes the zero-temperature phase transition in two dimensions. 

1. Renormalized classical regime of spin fluctuations. 

The spin susceptibility Xsp (q, 0) below Tx is almost singular at the antiferromagnetic wave vector Q2 = (tt, tt) 
because the energy scale 6U = Umf.c — Usp {Umf,c = 2/xo (Q, 0)) associated with the proximity to the SDW instability 
becomes exponentially small. |2^] This small energy scale, SU « T, leads to the so-called renormalized classical 
regime for the fluctuations. |]66| In this regime, the main contribution to the sum over Matsubara frequencies entering 
the local-moment sum rule Eq.(P8|) comes from iqn — and wave vectors (q — Q) < S,^^ near Q. Approximating 
Xsp (q, 0) by its asymptotic form 
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we obtain, in d dimensions 

where ct^ = n — 2(n|n|) — C < 1 is the left-hand side of Eq.(p8|) minus corrections C that come from the sum over 



2 2T /• A 1 ,„. 



non-zero Matsubara frequencies (quantum effects) and from (q — Q) ^ ^~^. There is an upper cutoff to the integral 
which is less than or of the order of the Brillouin zone size. The important point is that the left-hand side of the above 
equation Eq.(|5^) is bounded and weakly dependent on temperature. This implies, as discussed in detail in Ref. [ p3[ , 
that the above equation leads to critical exponents for the correlation length that are in the spherical model (n —^ 00) 
universality class. For our purposes, it suffices to notice that the integral converges even when ^ — > cx) in more than 
two dimensions. This leads to a finite transition temperature. In two dimensions, the transition temperature is pushed 
down to zero temperature and, doing the integral, one is left with a correlation length ^ that grows exponentially 
below Tx 

^^e^pl^na^eo^y (53) 

The important consequence of this is that, below Tx, the correlation length quickly becomes larger than the single- 
particle thermal de Broglie wave length S,th — vf/ (t^T). This has dramatic consequences on quasiparticles in two 
dimensions. 

2. Effect of critical spin fluctuations on quasiparticles. 

When the classical fluctuations {iqn — 0) become critical, they also give, in two dimensions, a dominant contribution 
to the self-energy at low frequency. To illustrate what we mean by the classical frequency contribution, neglect the 
contribution of charge fluctuations and single out the zero Matsubara frequency component from Eq. (E^) to obtain 
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Here, Ck is measured relative to the chemical potential. The last term is the contribution from quantum fluctuations. 
In this last term, the sum over Matsubara frequencies iqn must be done before the analytical continuation of ikn to 
real frequencies otherwise this analytical continuation would involve going through complex plane poles of the other 
terms entering the full sum over ig„. The contribution from classical fluctuations, iqn — 0, does not have this problem 
and furthermore it has the correct asymptotic behavior at ikn — > oo. Hence the contribution of classical fluctuations to 
the retarded self-energy S^ (k, lo) can obtained from the iqn = term by trivial analytical continuation ikn — > tj + iO. 
Note also that the chemical potential entering G'-"^ in the self-energy formula is /io = /^ = at half-filling. 

Doing the same substitution as above for the asymptotic form of the spin susceptibility Eq. (|5y) in the equation for 
the self-energy Eq.fifl) one obtains the following contribution to E from classical fluctuations 

Eez {K^kn) = ^J J^^TZ^^-i—j;^^^^-^:^^- (55) 

where we have expanded Ck-i-Q-i-q — Ck+q + q • Vk+q. In the case that wc consider, namely half- filling and k = hp, we 
have /io = /z = and ek^^+q = 0. The key point is again that in two dimensions the integral in this equation Eq.(|55|) 
is divergent at small q for ^ = 00. In a Fermi liquid, the imaginary part of the self-energy at the Fermi surface {uj = 0) 
behaves as S'j^(ki?, 0) ^ T^. Here instead, we find a singular contribution 

E'^(kf ,0) (X T [ d''-\^ , ] ex Te--" (56) 

J 91 + s 

that is proportional to ^ in d = 2 and hence is very large S^(ki?,0) « ~U^/{^th^'^) > 1 when the condition ^ > ^th 
is realized. By contrast, Ofor d = 3, E'^(ki;',0) ~ — t/(ln^)/ {^Q^th), so that the Fermi liquid is destroyed only in 
a very narrow temperature range close the Neel temperature T/v- Dimensional analysis again suffices to show that 
in four dimensions the classical critical fluctuations do not lead to any singular behavior. Three dimensions then is 
the upper critical dimension. As usual, logarithmic corrections exist at the upper critical dimension. The effect will 
be very small in three dimensions not only because it is logarithmic, but also because the fluctuation regime is very 
small, extending only in a narrow temperature range around the Neel temperature. By contrast, in two dimensions 
the effect extends all the way from the crossover temperature, Tx, which is of the order of the mean-fleld transition 
temperature, to zero temperature where the transition is. 

Wave vectors near Van Hove singularities are even more sensitive to classical thermal fluctuations. Indeed, near 
this point the expansion should be of the type eki/H+q+('r-Tr) ^ l'x~1y- This leads, in two dimensions, to even stronger 

divergence in S'^(ki?, 0) ex T^^ J dqy [(Sq^ + l) \qy\\ [Q. Even if the logarithmic divergence is cutoff the prefactor 
is larger by a factor of ^ compared with points far from the Van Hove singularities. 

B. Precursors of antiferromagnetic bands in two dimensions. 

Let us analyze in more details the consequences of this singular contribution of critical fluctuations to the self- 
energy in two dimensions. The integral appearing in the two-dimensional version of the expression for the self-energy, 
Eq. (p5|) , can be performed exactly |£j 

T.tKF,lkn) — I In „ h T< (57 

Here 7?. is a regular part. 

As a first application, we can use this expression to understand qualitatively both the temperature and size depen- 
dence of the Monte Carlo data for z{T) appearing in Fig. (2) of Ref. pQl or in the lower panel of Fig.(R) . Indeed, 
z{T) can be written as the alternating series —2G (ki?, /3/2) = — 4T^^^^ (—1)" / (^n — S" (k_Fj ikn))- Even though 
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the series converges slowly, in the beginning of the renormalized classical regime and for qualitative purposes it suf- 
fices to use the first term of this series. Then, using the expressions for the correlation length Eqs. ( p3| ) and for the 
self- energy (pj). One finds. 



z(T) 



rp2 



a^UUs: 



% , Tx-T^Tx. 



(58) 



On the infinite lattice, ^ starts growing exponentially below Tx, quickly becoming much larger than S,th- This implies 
z{T) ~ T^. On finite lattices f ~ ViV, which explains the size effect observed in Monte Carlo i.e. smaller z for smaller 
size N, {^th{Tx) - 5 for Fig.(|)). 

The analytic continuation of E {^f, ikn) in Eq.(p7[) is 
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For the wave vectors k away from the Fermi surface the anomalous contribution due to the classical fluctuation has 
a similar form but with lo replaced by (w — ek+q)- When T > Tx, the correlation length ^ becomes of order unity 
and, as we will show in Appendix (O), the regular part TZ dominates so that one recovers standard Fermi liquid 
behavior. Furthermore, even for large correlation length the regular part cannot be neglected when lo ^ T since the 
term exhibited here becomes small. Hence we concentrate on small frequencies and on T < Tx where the regular 
part TZ can be neglected. 

Exactly at the Fermi level [ui = 0) we recover the result of the previous section, namely that the imaginary part 
of the self-energy for ^ > £^th increases exponentially when the temperature decreases, E"(ki?,0) '-^ U£,/{^th£,o) oc 
T^ oc Texp (tto-^^qC/sp/T). The above analysis shows by contradiction that in the paramagnetic state below Tx there 
is no Fermi-liquid quasiparticle at kp, yet the symmetry of the system remains unbroken at any finite T. Indeed, 

starting from quasiparticles I Ga ) we found that as temperature decreases, E^(ki7, 0) increases indefinitely instead 

of decreasi ng.^ in direct contradiction with the starting hypothesis. By contrast, a self-consistent treatment where we 
use in Eq.([46|) the full Ga with a large S'j^(ki^, 0) shows that, for T < Tx , S'j^(ki?, 0) remains large in d = 2 an d doe s 
not vanish as T — > 0, again confirming that the system is not a Fermi liquid in this regime (See however Sec.(VIB) 
below). Strong modifications to the usual Fermi liquid picture also persist away from half-filling as long as Tx{n) > 0, 
as we discuss later. 

One can check that the large S]^(ki?,0) in two dimensions (for T < Tx) leads to a pseudogap in the infinite 
lattice, contrary to the conclusion reached in Ref. pq] . Indeed, instead of a quasiparticle peak, the spectral weight 
A(kF,i^) = — 2ImG_R (ki?, w) has a minimum at the Fermi level lo = and two symmetrically located maxima away 
from it. More specifically, ior vp/^ < \uj\ < T we have 
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The maxima are located aX lj — -ii^JUlIspa /2. These two maxima away from zero frequency correspond to precursors 
of the zero-temperature antiferromagnetic (or SDW) bands (shadow bands 133]). There is no quasiparticle peak b etwe en 
these two maxima when S,> ^th- This remains true in the case of no perfect nesting as well [g6[ (see also Sec.( VE )). 
We note that this is different from the results of the zero-temperature {£,th = oo) calculations of Kampf and Schrieffer 
p5[ that were based on a phenomenological susceptibility separable in momentum and frequency xk.sh. — fi'ijgi^)- 
As was explained in Ref. pq ] , the existence of precursors of antiferromagnetic bands (shadow bands in the terminology 
of Ref. [pq|) at zero temperature is an artifact of the separable form of the susceptibility. The third peak between the 
two precursors of antiferromagnetic bands that was found in Ref. [|35| is due to the fact that at zero temperature the 
imaginary part of the self-energy I]"(k, a; = 0, T = 0) is strictly zero at all k. In our calculations, precursor bands 
appear only at finite temperature when the system is moving towards a zero-temperature phase transition. In this 
case, the imaginary part of the self-energy goes to infinity for k on the "shadow Fermi surface" limr^o 5^"(k_F + Q, 0) oc 
T^ oc Texp (Gst/T) -^ oo and to zero at all other wave vectors. This is consistent with the SDW result which we 
should recover at T = 0. Indeed, the latter result can be described by the self-energy T,^{'k,u}) = A'^ / {uj — eCk+CD+ir]) 
which implies that the imaginary is a delta function S"(k,w) = — 7r^(a; — e(k-f Q)) instead of zero at all k as in a Fermi 
liquid. We note also that analyticity and the zero value of E" (k, a; = 0) in Ref. |35| automatically implies that the 
slope of the real part of the self-energy dT,'(k, uj)/duj\^^o is negative. By contrast, in our case dT^'ikp + Q, uj) / dujl^o^n 
is positive and increases with decreasing temperature, eventually diverging at the zero-temperature phase transition. 
The real part of the self-energy obtained using the asymptotic form Eq.dsS) is at the bottom left corner of Fig.(0) 
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with the corresponding spectral function A (kp, cu) shown above it. In Fig.(H) we have already shown the evolution of 
the spectral function A (kj?, a;) with temperature. The positions of the precursors of antiferromagnetic bands scale like 
a/2 which itself, at small coupling in two dimensions, scales like the mean field SDW transition temperature or gap 
(see Appendix B of Ref. |^^). As U increases, the predicted positions of the maxima obtained from the asymptotic 
form Eq. ( |60| ) will be less accurate since they will be at intermediate frequencies and the regular quantum contribution 
to the self-energy will affect more and more the position of the peaks. 

We have predicted l^^l that the exponential growth of the magnetic correlation length ^ below Tx will be accom- 
panied by the appearance of precursors of SDW bands in A(ki?,w) with no quasiparticle peak between them. By 
contrast with isotropic materials, in quasi-two-dimensional materials this effect should exist in a wide temperature 
range, from Tx (Tx «U < Ep) to the Neel temperature Tx [Tx - Tx ^ 10^ K). 

C. Contrast between magnetic precursor effects and Hubbard bands 

Although there are some formal similarities between the precursors of antiferromagnetic bands and the Hubbard 



bands (see Sec. Vl) we would like to stress that these are two different physical phenomena. A clear illustration of this 
is when a four peak structure exists in the spectral function ^(k, w), two peaks being precursors of antiferromagnetic 
bands, and two peaks being upper and lower Hubbard bands. The main differences between these bands are in 
the \i— dependence of the self-energy I](k, lj) and in the conditions for which these bands develop. Precursors of 
antiferromagnetic bands appear even for small U in the renormalized classical regime T < Tx , and their dispersion has 
the quasi-periodicity of the magnetic Brillouin zone. In contrast, upper and lower Hubbard bands are high-frequency 
features that appear only for sufficiently large U > W and T < U and have the periodicity of the whole Brillouin 
zone in the paramagnetic state. Furthermore, the existence of Hubbard bands is not sensitive to dimensionality so 
they exist even in infinite dimension where the self-energy does not depends on moment um k at all. In contrast, the 



upper critical dimension for the precursors of antiferromagnetic bands is three (see Sec.(Vn)). 



In our theory the precursors of antiferromagnetic bands come from the almost singular behavior of the zero Matsub- 
ara frequency susceptibility Xsp (q, 0), which leads to the characteristic behavior of S(k, w) = A|^ ^/ (w — £(k -I- Q)) 
with A|,j g ex Tln(^). On another hand, the Hubbard bands appear in our theory because the high-frequency asymp- 
totics S(k, (jj) ex Afj g/uj has already set in for uj > W, and this leads to the bands at cj = ±Ah.b for A > W. 
(see for more details Sec.(|Vl|)). The coefficient A^ g is determined by the sum over all Matsubara frequencies and q: 
Ajjg = TUN-^ X^q,™ [C^spXsp(q, i^n) + UchXchiq., i^n)]- 

It was noticed in Monte Carlo simulations [Q [p8| that for intermediate C/, the spectral weight has four maxima. We 
think that peaks at a; ~ ±[//2 are Hubbard bands, while the peaks closer to w = are precursors of antiferromagnetic 
bands. If this interpretation is correct, then the latter peaks should disappear with increasing temperature when ^ 
becomes smaller than ^th, while the Hubbard bands should exist as long a,s T < U. 

While the location of the precursors of antiferromagnetic bands should be accurate in our theory, the same will 
not be true for the location of the upper and lower Hubbard bands. This is because our theory is tuned to the low 
frequency behavior of the irreducible vertices and does not have the right numerical coefficient in the high-frequency 



expansion of the self-energy, as shown in Eq.(ElC) below. Nevertheless, our analytical approach to date is the only 
one that agrees at least qualitatively with the finding that precursors of antiferromagnetic bands as well as upper and 
lower Hubbard bands can occur simultaneously. Note however that a four peak structure at n = 1 was also obtained 
in Ref. iTfl] but the physical difference between Hubbard bands and precursors of antiferromagnetic bands was not 



clearly spelled out. We comment on recent findings of the FLEX approach in Sec.(V]) [p9| |3 



D. Can the precursors of antiferromagnetic bands exist in three dimensions? 

In two dimensions, the finite-temperature phase is disordered, but the zero-temperature one is ordered and has 
a finite gap, except at the quantum critical point away from half-filling. Hence, precursors of antiferromagnetic 
bands that appear in the paramagnetic state do so with a finite pseudogap which appears consistent with the finite 
zero-temperature gap towards which the system is evolving. By contrast, in higher dimensions the gap opens-up 
with a zero value at the transition temperature. Based on this simple argument, one does not expect precursors of 
antiferromagnetic bands in dimensions larger than two (see, however, below). Here, we will also show that there is no 
phase space reasons for the existence of precursors of the antiferromagnetic bands when d > 2. 

We have already shown that in three dimensions the quasiparticle at the Fermi level at half-filling will have an 
imaginary part of the self-energy that grows like Tln^, an effect that is much weaker than T^ found in two dimensions. 
Despite this small effect, in three dimensions the classical fluctuations do not affect the self-energy for energies larger 
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than vpS^ ^. Indeed, consider the contribution of classical thermal fluctuations to the self-energy Eg. (|55|) . In two 
dimensions, we have for \uj\ > vp^~^ 

ReRf(k„„,].g/||,-^i. ,61) 

which allows us to recover the approximate formula for the spectral weight given in Eq.(pG) above. In three dimensions 
however, this approximation cannot be done because the integral is not dominated by small values of q. To see this 
explicitly in three dimensions, consider the contribution of classical thermal fluctuations 
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As long as |w| > vf£,~^, the logarithmic singularity that develops at q\\ = when ^^^ -^ is integrable and gives no 
singular contribution to the self-energy. Hence, unusual effects of classical thermal fluctuations are confined to the 
range of frequencies |a;| < vp£,~'^ . At higher frequencies, |w| > wfC~^i all bosonic Matsubara frequencies in Eq.(El) need 
to be taken into account and from phase space considerations alone there is no reason for the existence of precursors 
of antiferromagnetic bands in the iD case. However, the existence of such bands in 313 cannot be completely excluded 
based on dimensional arguments alone because they occur at finite frequencies and strictly speaking they are non- 
universal. In particular, as discussed in Ref. p3], one expects to see precursors that look like 2D antiferromagnetic 
bands (shadow bands) in the vicinity of the finite temperature phase transition in strongly anisotropic quasi-two- 
dimensional material. On the other hand, such bands do not generically exist in the almost isotropic 3-D case, because 
even in 2D the conditions for such bands are quite stringent. The difference between shadow bands and Hubbard 
bands has been discussed in the previous subsection and the discussion of non-analyticities sometimes encountered in 
Fermi liquid theory can be found in Appendix (O). 

E. Away from half-filling 

Close to half-filling, in the nearest-neighbor hopping model, one can enter a renormalized classical regime with 
large antiferromagnetic correlation length, even though the zero-temperature Fermi surface properties may favor 
incommensurate correlations. This renormalized-classical regime with large (tt, tt) correlations occurs when Tx 3> /io- 
By arguments similar to those above, one finds that in this regime one still has precursors of antiferromagnetic 
bands. However, the chemical potential is in or near the lower precursor band and the system remains metallic. The 
high-frequency precursor appears only below Tx at w w Ek+q- 

With second-neighbor hopping, the points of the Fermi surface that intersect the magnetic Brillouin zone (hot 
spots) behave as does the whole Fermi surface of the nearest-neighbor (nested) case discussed above. These questions 
were discussed in detail in Ref. |p6|| . 

F. The pairing pseudogap and precursors of superconducting bands in two dimensions 

As we have already pointed out above, the results for the single particle spectra obtained for the half-filled nearest- 
neighbor hopping repulsive Hubbard model can be directly applied to the corresponding attractive Hubbard model, 
in which case the pseudogap opens up in the renormalized classical regime of pairing and charge fluctuations. Away 
from half-filling, the symmetry between charge and pair correlations is lost and pair fluctuations dominate, becoming 
infinite at the Kosterlitz-Thouless transition temperature. This temperature is below the temperature at which the 
magnitude of the pair order parameter acquires rigidity despite the randomness of its phase. One expects then that 
a pseudogap will also open in this case when the correlation length for pairing fluctuations becomes larger than 
the single-particle thermal de Broglie wavelength ^pairing > £,th = vf/T. This should occur below the crossover 
temperature to the renormalized classical regime of pairing fluctuations but above the Kosterlitz-Thouless transition 
temperature. 

The quantitative microscopic theory for the negative U Hubbard model will be considered in a separate publication. 
By contrast with all other sections of this paper, our considerations here will be more phenomenological. Nevertheless, 
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they will allow us to present some analytical results for the self-energy obtained in the critical regime dominated by 
pairing fluctuations. Details of the model should not be very important since we are in a regime where everything is 
dominated by long wave length fluctuations. 

The derivation of I](k , m) in t he pairing case is a straightforward extension of what we did in the antiferromagnetic 
case (see Subsecs.(VA2), (VB) and Ref. [p6|). In particular, in complete analogy with the magnetic case, the main 
contribution to the self-energy in the critical regime comes from the classical thermal fluctuations iqn — 0. Assuming 
some effective coupling constant g' between quasiparticles and pairing fluctuations, which in general can be momentum 
dependent, one can write in the one loop approximation 
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Here et is the electron dispersion relative to the chemical potential, and all factors in front of integral are reabsorbed 
into the coupling constant g' . This expression is similar to the expression Eq.(|5^) in the magnetic case but there 
are two important differences: i) instead of ek+Q+q we have now e_k+q; n) there is no minus sign in front of e_k+q- 
The first difference is due to the fact that superconductivity usually occurs with zero center of mass momentum for 
the pair, and hence the pairing susceptibility in the normal state Xp °c ^/i^.p'^ + Q^) must be peaked near q = 0, 
(the integration variable q in Eq. ( ^5|) was measured relative to Q = (tt, tt)). The second difference comes from 
the fact that we are now considering the contribution to E coming from the particle-particle channel instead of the 
particle-hole channel. Taking the integrals over q and using the fact that small q only will contribute we neglect the 
q dependence of the coupling constant and obtain for the imaginary part of T.ci the following expression 
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In the renormalized classical regime the pairing correlation length ^p increases faster with decreasing temperature 
than ^th = vf/T. Consequently, S"(ki?,0) tends to diverge with decreasing temperature and a pairing pseudogap 
in the spectral function A{\i.p^uj) opens up over the complete Fermi surface, except maybe at a few points where 
.g'(k) = 0. This is different from the antiferromagnetic case, where the pseudogap in ACkp = ^h.sp.,'-^) opens up 
only when, so called, "hot spots" {eCkh.sp. + Q) = ^(kh.sp.) = 0) exist in a given model Bq]. The antiferromagnetic 
pseudogap opens everywhere on the Fermi surface only in the case of perfect nesting, where all points on the Fermi 
surface are "hot spots" . 

The real part of the self-energy can be obtained from Eq. (pq) using the Kramers-Kronig relation and has the form: 
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To understand how precursors of the superconducting bands develop, let us look at S'(k, u) at frequencies 
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-k^p • In this case, using inversion symmetry £_k — £k, one can obtain from Eq.(66) the following asymptotic form 
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When, ^p r^ exp{const/T) (see, more general case below) this form of the self-energy leads to the usual BCS result 
T.'{k,uj) w A2(k)/(w-hek) with the gap A^ (k) w (g'(k)/27r)Tln^p. On the other hand, the imaginary part S"(k,a;), 
Eq. (pSi) , vanishes everywhere in the T = limit, except when w = —i^ where it becomes infinite. The results for S' 
and E" can thus be combined to write for the corresponding limit of the retarded self-energy E^ = A^ (k) /{Lu+e\^+ir]). 
This limit leads to the standard BCS expression for the normal Green's function when substituted back into the Dyson 
equation G^ = 1/ (cj -I- «?? — Ck ~ E^(k, a;)). Above the transition temperature, the anomalous Green's function 
remains zero since there is no broken symmetry. The qualitative picture for the development of the pairing pseudogap 
and of the precursors of superconducting bands at k = k^ is illustrated in Fig.(g) and in the left part of Fig.([ll|). 
While in the case of magnetic critical fluctuations these figures describe the precursor effect in ACkp,^) for perfect- 
nesting or for the "hot spots" (when such points exist), in the case of pairing fluctuations they describe the spectra 
for all hp and for all fillings where the ground state is superconducting. 

We need to comment on a subtle difference between the antiferromagnetic and the pairing precursor effects in 
the single particle spectra. While the magnetic order parameter has three components and can order only at zero 
temperature in the two-dimensional repulsive model, away from half-filling in the attractive model the pairing order 
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parameter becomes the only relevant order parameter at low temperature. Since it has only two components, a finite 
temperature Kosterlitz-Thouless phase transition is then allowed in two dimensions. The critical behavior in vicinity 
of this transition is given by ^p oc ex])[const/{T — TktY^^] instead of ^ ex exp{const/T) as in the magnetic case. 
To take this properly into account one would need a treatment of the problem that is more sophisticated than that 
given above. In particular, one would have to take into account corrections to the simple form that we used for 
the pairing susceptibility Xpi.<lj 0) oa l/(^j7^ + <Z^)- This Lorentzian form of the susceptibility in the critical regime is 
strictly valid only in the n — oo limit ( n is the number of the components of the order parameter) and is, clearly, 
a less accurate approximation in the case of pairing fluctuations {n — 2) than in the case of the antiferromagnetic 
fluctuations {n ~ 3). Nevertheless, we believe that qualitatively the picture given above is correct for two reasons. 
First, because in the Kosterlitz-Thouless picture the magnitude of the order parameter is locally non-zero starting 
below a crossover temperature Tx that is larger than the transition temperature Tr-t. It is only the phase that 
is globally decorrelated above Tkt- This means that locally the quasiparticles are basically in a superconducting 
state even above Tkt- A second reason to believe in the precursor effects is that the superfluid density and the gap 
are finite as T — j Tj^ j, and, hence, the two peak structure in A(ki^,w) exists even as the phase transition point 
is approached from the low-temperature side. By analogy with the antiferromagnetic case, this two peak structure 
should not immediately disappear when one increases the temperature slightly above Tkt- 

Finally, we point out that the precursor phenomenon described above has to be distinguished from, so-called, pre- 
formed pairs considered first by P. Nozieres and S. Schmitt-Rink |7l| (see also |Q). These pre- formed pairs exist 
in any dimension when the coupling strength is sufficiently large, while the precursor effect considered above can be 
caused by arbitrarily small attractive interactions but only in two dimensions. We think that recent Monte Carlo 
data [Q on the negative U = —W/2 Hubbard model illustrates the opening of the single-particle pseudogap due to 
critical fluctuations, rather than a strong-coupling effect. In these simulations, the drop in the density of states at 
the Fermi level should be accompanied by a simultaneous rapid increase of the pairing structure factor S'p(q =0, T) . 
The latter must be exponential in the inflnite 2D lattice and a size analysis of Monte Carlo data similar to the one 
shown on Fig.{0) would be extremely helpful to clarify this issue. 

VI. ABSENCE OF THE PRECURSORS OF ANTIFERROMAGNETIC BANDS AND UPPER AND 
LOWER HUBBARD BANDS IN ELIASHBERG-TYPE SELF-CONSISTENT THEORIES 

In this section, we explain why the theories that use self-consistent propagators but neglect the corresponding 
frequency-dependent vertex corrections fail to see two important physical effects: namely upper and lower Hubbard 
bands, as well as the precursors of antiferromagnetic bands that we just discussed. The failure of this type of self- 
consistent schemes to correctly predict upper and lower Hubbard bands has been realized a long time ago in the 



context of calculations in inflnite dimension |74[| 1 11 1 . While one may brush aside this failure by claiming that high- 
energy phenomena are not so relevant to low-energy physics, we show that in fact these schemes also fail to reproduce 
the low- energy pseudogap and the precursors of antiferromagnetic bands for essentially the same reasons that they 
fail to see Hubbard bands. It is thus useful to start by a discussion of the better understood phenomenon of upper 
and lower Hubbard bands and then to move to precursors of antiferromagnetic bands. 

A. Why Eliashberg-type self-consistency for the electronic self-energy kills Hubbard bands 

We first note that ordinary perturbation theory satisfies the correct high-frequency behavior Eq.(p^) for the self- 
energy namely, for fc„ ^ W 

hm S^ (k,zfc„) ^ [/n_. + ^'"'" ^} ~ "'"^ + . . . (68) 

It is the latter property that guarantees the existence of the Hubbard bands for U > W. To see this, consider the 
half-filled case. In this case, n^a- — 1/2, /i = U/2 and one finds for the spectral weight 



l^ 4c. 
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Ai^>-)-Z ^77^—- (69) 



which has pronounced maxima at the upper and lower Hubbard bands, namely u — ±t//2,has long as S"is not too 
large. Since these results are obtained using high-frequency asymptotics, they are valid only when the asymptotic 
Eq.(|68[) has already set in when uj ^ C//2.In the exact theory and in ordinary perturbation theory in terms of bare 
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Green functions G^^\ Eq.(|6^) is valid for \uj\ ^ M^and the Hubbard bands appears as soon as t/become larger than 
W. 

The fact that this simple high-frequency behavior sets in at the energy scale given by M^rather than U, even when 
W < U, is a non-trivial consequence of the Pauli principle. To see this we first recall the exact result for the self-energy 
So- (k,ifc„)in the atomic limit H 

E^— (k,zfc„) = t/n_o + /^'"-'^(l-"-'^) (70) 

tkn + n - U [1 - n-a) 

Formally, the atomic limit means that hopping is the smallest of all energy scales in the problem, including the 
temperature, t <^ T, which is not a very interesting case. However, the same arguments that have been used to 
derive the expression ( [70| ) in the atomic limit can be used to show that Eq. ([70|) is valid at any T/iwhen fc„ ^ W. 
Indeed, in the equations of motion for two-particle correlators W] one can neglect hopping terms when A:„ ^ W. This 
is where the asymptotic behavior (|70|) sets in since the equations of motion then immediately lend themselves to a 
solution without any additional approximation for the interacting term. This solution is possible because the Pauli 
principle nf^ — TT-io-allows us to collapse three-particle correlation function which enters equation of motion to the two- 
particle one U (Tr (rii-a {n) iH-a (n) Cia (n) c]o (■^i) ) / = U (Tt (rii-a (n) Cia {n) c]^ (■^i) ) /■ Hence, the expression 
for atomic limit Eq.(|7^) is also a general result for the self-energy that is valid for fc„ ^ W.At half-filling n_o- — 1/2, 
/i — t//2and the asymptotic ^68) sets in at fc„ ~ W, as was pointed out above. Away from half-filling, as long 
as l/i — S (cxD)|and \fi — U {1 — n_CT)|are both much smaller than T/F,(they both vanish at half-filling), the asymptotic 
behavior will also start at fc„ ~ W. 

The situation is qualitatively different when one uses dressed G reen functions, but does not take into account the 



frequency dependence of the vertex, as it is done in FLEX (See Eq. E9)) or for second-order perturbation theory with 
dressed G. For example, the second-order expression for E^- (k,z/c„)in terms of full Gdoes satisfy the asymptotics 
Eq. (|68|) , but it sets in too late, namely for fc„ » U, instead of fc„ :^ W. Indeed, when fc„ ^ W, the equation for the 
self-energy at half-filling in this type of theories reduces to 

where A^ ~ cC/^/4with ca constant of proportionality involving the sum over all wave vectors and Matsubara fre- 
quencies of the self-consistent dynamical susceptibilities. In a given theory the value of c may differ from its value 
c = lobtaincd from the exact result Eq. d7Q), but its always of order unity. The solution of Eq.(|7l|) 



E {ikn) = \ikn - i \j{iknf - 4A2 (72) 

has the analytically continued form 

ReE« (cc) = I - ^^ (kl - 2A) Vc.2_4A2 (73) 

ImE^^ (w) = --e (2A - \uj\) V4A2 - uj^ (74) 

From this one can immediately see that a f/^/tjregime exists for ReE-'^ (ti;)only when \uj\ ^ [/,(with 2A = U). 

This means that such regime sets in too late to give the Hubbard bands described by Eq. ( |69| ) , because the Hubbard 
bands occur at w = ±C//2 and for such o^the asymptotic form E^ oc U'^ /ujis not valid yet in FLEX and similar 
theories. Consequently, instead of well defined peaks at a; = ±C//2 in the half-filled case, one obtains only long tails 
in the spectral function A^, (k,a;) , no matter how large U is. [Q(see also following subsection). 

This explains why there is no Hubbard bands in any theory that uses self-consistent Green's functions, but neglects 
the frequency dependence of the vertex. This is an explicit example that illustrates what seems to be a more general 
phenomenon when there is no Migdal theorem for vertex corrections: a calculation with dressed Green's functions 
but no frequency dependent vertex correction often gives worse results than a calculation done with bare Green's 
functions and a frequency independent vertex. 
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B. Why FLEX fails to see precursors of antiferromagnetic bands 

In this subsection we describe the quaUtative differences between our resuhs and the results of FLEX approximations 
given by Eq. (^) with regards to the "shadow bands" and explain why we believe that the failure of the FLEX to 
reproduce these bands is an artifact of that approximation. To avoid any confusion, we first clarify the terminology, 
because the term "shadow bands" has been used previously to describe different physical effects, (see for details Ref. 



|36|). We note that the so-called shadow features discussed in |36|, ||33] as well as the pseudogap in the total density 
of states N{u!) = {1/N) J^w -^ 0^: i^))exist in both theories and we will not discuss them here. Instead, we concentrate 
on the precursors of antiferromagnetic bands in the spectral function A (kj?, lj) which correspond to two new solutions 
of the quasi-particle equation 

w-e(k)+/x-E(k,cj) ==0 (75) 

We start by recalling a simple physical argument why the precursors of antiferromagnetic bands must exist at finite 
temperatures in the vicinity of the zero-temperature phase transition in two dimensions. This can be best understood 
by contrasting this case with isotr opic 3 D case where such precursor effect arc highly unlikely. (For a discussion of the 



strongly anisotropic case see Sec.(VD)). Indeed, in three dimensions there is a finite temperature phase transition 
and the gap is equal to zero at this temperature A(TAr) = 0. Consequently at T^rthere is only one peak in the 
A{]<ip,u))a,i to = Owhich starts to split into two peaks only below Tjv. Based on this simple physical picture, one 
would not expect to see precursors of antiferromagnetic bands above T^rin this case. The situation is qualitatively 
different in two-dimensions where classical thermal fluctuations suppress long-range order at any finite temperature 
while at the T — phase transition the system goes directly into the ordered state with a finite gap. Clearly, the two 
peak structure in A (\<.p,Lu) at T = cannot disappear as soon as we raise the tem perature. 



For simplicity we again consider half- filling. As we have seen in section (VB) two new quasi-particle peaks do 
appear in the renormalized classical regime T < Txin our theory. We have also found a pseudogap with the minimum 
at w = Oin this regime. In contrast, the numerical solution of the FLEX equations |3q] found a spectral function with 
a single maximum in A{'kp,ijj)a.t to = Oeven when x-'^"'^'^(q, O)becomes strongly peaked at q = Q. With decreasing 
temperature this central maximum becomes anomalously broad, but the two peak structure does not appear. The 
clear deviation from the Fermi liquid is signaled by the positive sign of 9E' {\<ip,ijj) jdio > 0. However the value of 
dT,' (kp,uj) /diudocs not become larger than unity. The latter would unavoidably lead to the existence of two new 
quasi-particle peaks away from w = Oas is clear from the graphical solution of the quasiparticle equation Eq.(|75|) 
shown on the bottom left panel of Figs.(pT|). 

We now explain analytically the origin of these qualitative differences in the two theories. In our theory 
dY}' {hp , uj) / dLLi\uj=o oc T^^and hence it quickly becomes larger than unity in the renormalized classical regime 
^ ex e'xp{const/T). In addition, for lj > vp£^~^the real part of the self-energy has the same behavior as in the 
ordered state E(k^,a;) oc A^/wwith A^ c>c Tln^ — const. The important point is that this asymptotic behavior 
T,{'kp,Lu) (X A^/cjof the self-energy already sets in for cu ^ A ^ vp^~^. It is this property that leads to the ap- 
pearance of the precursors of antiferromagnetic bands at lo — ±Ain a manner analogous to the appearance of the 
Hubbard bands in the strong coupling limit that is discussed in the previous subsection. Let's now try to understand 
analytically what happens in the FLEX approximation. As in our theory, the main contribution to the self-energy 
in the strongly fiuctuating regime comes from the zero-frequency term in the Matsubara sum in the equation for the 
self-energy (Eq.(|54|) in our theory and Eq.dE9|) in FLEX). An upper bound of the effect of the critical spin fluctuations 
can be obtained by approximating Tx^^'^(q, 0) ex ^(q). Then one immediately obtains the same expression for the 
self-energy as the one obtained in FLEX in the context of Hubbard bands Eq.([7l|). (The only difference is that the 
parameter A is now defined by the zero- frequency Matsubara contribution of x^^^(flj 0), rather than by the sum over 
all Matsubara frequencies.) As we have already discussed in the context of Hubbard bands, such a form for S does 
not lead to the appearance of two new quasiparticle solution away from uj = because the characteristic behavior 
E(ki^,(jj) ex A'^/lu sets in too late, namely for w ^ A. In addition, the slope of E' (kp,uj) at u; ^ does not diverge 
with decreasing temperature as in our theory but instead saturates to its value given by the analog of Eq. ([73|) , i.e. 
9E' (kpjLLi) /du) < 1/2. As we mentioned above, a value larger than unity 9E' (kpjU!) jdio > 1 would guarantee the 



existence of two new solutions of the quasiparticle equation Eq. ([75[) away from uj = 0. The right-hand side of Fig.(|ll|) 



illustrates clearly what happens in a FLEX-like approach such as Eq.(71). The contribution of classical fluctuations to 
the spectral weight does not lead to a Fermi liquid since A {kp,io) saturates to a finite width as temperature decreases, 
but nevertheless precursors of shadow bands do not occur because dT,' (k^, w) /duj is bounded below unity. (Note 
that the spectral weight would not vanish so steeply at large frequencies if we had taken into account the quantum 
contribution of the spin fluctuations, as in full FLEX calculations.) 

We just saw that the self-consistency in the propagators without corresponding self-consistency in the vertices 
inhibits the existence of the shadow bands in essentially the same way as it inhibits the existence of the Hubbard 
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FIG. 11. Top two panels are qualitative sketches of the spectral weight at the Fermi wave vector at half-filhng. The plots 
are obtained by retainig only the classical contribution to the self-energy for T < Tx using parameters corresponding to the 
typical (7 = 4, of Monte Carlo simulations. The two bottom panels are the corresponding plots of ReT,{uo). The left-hand side 
of this figure is obtained using our approximation while the right-hand side is obtained from the FLEX-like approach. The 
intersection with the 45 degree line uj in the bottom-left panel gives rise to the precursors of antiferromagnetic bands seen right 
above it. 



bands. It thus seems to us very likely that the absence of the precursors of antiferromagnetic bands below Tx in FLEX 
is an artifact. This conclusion can be reliably verified by comparison with Monte Carlo data despite the f act th at the 
latter is done for finite lattices and in the Matsubara formalism. This was discussed in more detail in Sec.( |VBD . Here 
we just note that the temperature dependence of Matsubara quantities such as G(ki?, t = (3/2) and E(ki?, iki) have 
a very characteristic form in the pseudogap regime. For example, S(k, iki) oc l/(i7rT) in the pseudogap regime, while 
in FLEX we would expect a much weaker temperature dependence of this quantity (the upper bound being given by 
the analog of Eq.Q). 

We also would like to comment on the ID model |75| ] which describes the interaction of electrons with static spin 
fluctuations characterized by the susceptibility Xsp oc ^(aj)[^~^/(q — Q)^ + £,^^]- The nice thing about this model is 
that it has an exact solution which shows the development of shadow bands and of the pseudogap in ACkpjOj). A 
treatment similar to ours which uses non-interacting Green's functions in the one-loop approximation also reproduces 
this feature JT^]. However, the analogous approximation with dressed Green's functions leads to Eq.(^l|) and hence 
inhibits the existence of the "shadow bands" and of the pseudogap in A (k^?, w). 

In closing we comment on semantics and on the physical interpretation of some results obtained in the FLEX 
approximation. The expression "conserving approximation" has been widely used to describe FLEX calculations of 
the single particle properties and, in particular, in the context of the shadow bands and of the failure of Luttinger's 
theorem [ p8|j3^j69|] . The conserving aspect has been emphasized, but in fact the only desirable feature in the calculation 
of the single-particle properties is that the self-energy E is obtained from a functional derivative of the Luttinger- 
Ward functional E = d^/SG and hence it is guaranteed to satisfy Luttinger's theorem whenever appropriate. Only 
on the next level does this scheme lead to a calculation of the "true" susceptibilities p4| and of collective modes 
that satisfy conservation laws (Ward identities). However, these "true" susceptibilities are never substituted back 
in the calculation of the self-energy and the effect of "true" collective modes on the single-particle spectrum is an 
open question in FLEX. In fact, the RPA propagators xbpa appearing in the self-energy expression are different 
from susceptibilities from which collective modes should be computed and further they explicitly break conservation 
laws, as can be seen from the fact that RPA-like exp ressio ns xr pa — Xo/(l — Uxo) with a dressed bubble xo have 
the unphysical properties that are mentioned in Eqs.( A2^ ) and (A27) of Appendix (Q). The fact that there are in 
effect two susceptibilities in the FLEX approximation leads, in our opinion, to some confusion and incorrect physical 
interpretation of the results in the literature. In particular, it was argued that the non-Fermi-liquid behavior and 
deviations from Luttinger theorem found in FLEX [ p8| , p7|j69| ] are not due to critical thermal fluctuation in the vicinity 
of the phase transition but are rather the result of large U. The reasoning for such claim was that although the RPA 
susceptibilities xrpa is very strongly peaked at q = Q, the "true" FLEX susceptibility is not. In our opinion, such 
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claim could be justified only if one would substitute the "true" susceptibility back in the calculation of E (for example 



using the exact Eq.(31)) and found that the deviation from the Luttinger theorem and other qualitative changes 
in A{k,LL!) increase with decreasing temperature without almost divergent behavior of the conserving susceptibility 
XspiQ, 0) and of the static structure factor Sgp (Q) ■ 

The Monte Carlo data on Fig. (M) are also instructive since they clearly show that qualitative changes in the single- 
particle spectra occur when the system enters the renormalized classical regime with rapidly growing Ssp (Q). The 
fact that the FLEX "true" susceptibility does not show such behavior at half-filling pSJ tells us that it even more 
drastically disagrees with the Monte Carlo data than the RPA-like x which enters the expression for self-energy. 
Moreover, even away from half-filling the "true" susceptibility in FLEX at q = Q significantly underestimates the 
strength of the spin fluctuations, as is clear from the comparison with Monte Carlo data in Fig.(g). In our opinion the, 
so-called, "true susceptibility" in FLEX is the key element in the confusion surrounding the interpretation of FLEX 
results for the self-energy because the "true susceptibility" never comes in the calculation of the self-energy For all 
practical purposes these calculations of the self-energy should be considered as consistent with Luttinger's theorem 
at T = but based on a non-conserving susceptibility. Consistency with conservation laws and consistency with 
Luttinger's theorem are not identical requirements because to satisfy rigorously Luttinger's theorem one needs that 
E = 6^/dG, while to have conserving susceptibilities one needs that the irreducible vertices used in Bethe-Salpeter 



equation Eq.(|2^) should be obtained from F = S'^^/SGSG. 

VII. DOMAIN OF VALIDITY OF OUR APPROACH 

Our approach is not valid beyond intermediate coupling. That is perhaps best illustrated by Fig.(y) that shows that 
the crossover temperature first increases with U and then saturates instead of decreasing. The decrease is expected 
on general grounds from the fact that at strong coupling the tendency to antiferromagnetism should decrease roughly 
as J ~ t^ /U. The reason for this failure of our approach is clear. As we know from studies in infinite dimension pi[ , 
to account for strong-coupling effects it is necessary to include at least a frequency dependence to the self-energy and 
to the corresponding irreducible vertices. 

Our theory also fails at half-filling deep in the renormalized classical regime, i.e. T <C Tx mainly for two reasons. 
First, the ansatz Ugp — [/^^^(O), Eq.(^), fails in the sense that 5|j,(0) eventually reaches zero at T — > because of the 
log T divergence in the irreducible susceptibility xo (tJ", tt) due to perfect nesting. The physically appropriate choice 
for (7tx(0) in the renormalized classical regime is to keep its value fixed to its crossover-temperature value (See Fig.(H) 
and Sec.([[V|)). The more serious reason why our approach fails for T <^ Tx is that, as we just saw, critical fluctuations 
destroy completely the Fermi liquid quasiparticles and lead to a pseudogap. This invalidates our starting point. It 
is likely that in a more self-consistent theory, the logarithmic divergence of the appropriate irreducible susceptibility 
will be cutoff by the pseudogap. However, just a simple dressing of the Green's function is n ot the correct solution to 



the problem because it would make the theory non-conserving, as we discussed in Sec. (A 3). One needs to take into 
account wave vector and energy dependent vertex corrections similar to those discussed by Schrieffer | f76[ [ [77| . 

VIII. COMPARISONS WITH OTHER APPROACHES 

In Appendix (W), we di scus s in detail various theories, pointing out limitations and advantages based on the criteria 
established in Appendix ( |A 2[ ) and ([A3| ). More specifically, we include in our list of desirable properties, the local 

Paifli principle («■?) = (^t) ' ^^^ Mermin- Wagner theorem Eq.( [Al7[ ), the Ward identities Eq.(A31), and /—sum rule 



Eq.( A25 ), one-particle vs two-particle consistency Eg. (l, l) Ga (l, 1^) = U {n-\ni) Eq. ([44|) , Luttinger's theorem, and 



the large frequency asymptotic for the self-energy Eq.(|68D, which is important fo r th e existence of the Hubbard bands. 



Comments on the forward scattering sum-rule Eq.([A3q) appear in Appendix (A 4). In the present section, we only 
state without proof where each theory has strengths and weaknesses. 

In standard paramagnon theories, pq ] [ p2[ the spin and charge fiuctuations are computed by RPA, using either bare 
or dressed Green's functions. Then the fluctuations are feedback in the self-energy When RPA with bare Green's 
functions are used for the collective modes, these satisfy the /—sum rule, but that is the only one of our requirements 
that is satisfied by such theories. 

In conserving approximation schemes p6| , p4| the Mermin- Wagner theorem, the Luttinger theorem and conserva- 
tion laws are satisfied, but none of the other above requirements are fulfilled. 

In the parquet approach, p5| | p3| one enforces complete antisymmetry of the four point function by writing down 
fully crossing-symmetric equations for these. However, in actual calculations, the local Pauli principle, the Mermin 
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Wagner theorem, and the consistency between one and two particle properties are only approximately satisfied, while 
nothing enforces the other requirements. 

In our approach, the high-frequency asymptotics and Luttinger's theorem are satisfied to a very good degree of 
approximation while all other properties in our list are exactly enforced. Let us specify the level of approximation. 
Luttinger's theorem is trivially satisfied with our initial approximation for the self-energy Eo- , but at the next level 



.(1) 



.(1) 



fixed. With 



of approximation, S^ , one needs a new chemical potential to keep the electron density Tr Ga (1,1) 

this new chemical potential the Fermi surface volume is preserved to a very high accuracy. Finally, consider the high- 
frequency asymptotics. Since we use bare propagators, the high-frequency asymptotics comes in at the appropriate 
frequency scale, namely ifc„ '^ W , which is crucial for the existence of the Hubbard bands. However, the coefficient of 
the 1/ikn term in the high-frequency expansion Eq.(pq) is incorrect because our irreducible vertices Usp and Uch are 
tuned to the low frequencies. If one would take into account the frequency dependence of Usp and Uch and assume 
that at high frequency they become equal to the bare interaction U, then one would recover the exact result, provided 
the Pauli principle in the form of Eq. (JS^) is satisfied. The difficulty with such a procedure is that frequency dependent 
irreducible vertices requires frequency dependent self-energy in the calculation of collective modes and that would 
make the theory much more complicated. Yet it is, probably, the only way to extend the theory to strong coupling. 

IX. CONCLUSION 

We have presented a new simple approach [g9[ jsy] to the repulsive single-band Hubbard model. We have also 
critically compared competing approaches, such as paramagnon, fluctuation exchange approximation, and pseudo- 
potential parquet approaches. Our approach is applicable for arbitrary band structure [ p4[ and gives us not only 
a quantitative description of the Hubbard model, but also provide us with some qualitatively new results. Let us 
summarize our theory again. We first obtain spin and charge fluctuations by a self-consistent parameterization of 
the two-particle effective interactions (irreducible vertices) that satisfies a number of exact constraints usually not 
fulfilled by standard diagrammatic approaches to the many-body problem. Then the influence of collective modes 
on single-particle properties is taken into account in such a way that single-particle properties are consistent with 
two-particle correlators, which describe these collective modes. More specifically, our approach satisfies the following 
constraints: 

1. Spin and charge susceptibilities, through the fiuctuation-dissipation theorem, satisfy the Pauli principle in the 
form (ni) — (n|) as well as the local moment sum-rule, conservations laws and consistency with the equations 
of motion in a local-field-like approximation. 

2. In two dimensions, the spin fluctuations satisfy the Mermin- Wagner theorem. 

3. The effect of collective modes on single-particle properties is obtained by a paramagnon-like formula that is 
consistent with the two-particle properties in the sense that the potential energy obtained from Tr [EG] is 
identical to that obtained from applying the fluctuation-dissipation theorem to spin and charge susceptibilities. 

4. Vertex corrections are included not only in spin and charge susceptibilities {Usp ^ Uch 7^ U) but also in the self- 
energy formula. In the latter case, this takes into account the fact that there is no Migdal theorem controlling 
the effect of spin and charge fluctuations on the self-energy 

The results for both single-particle and two-particle properties are in quantitative agreement with Monte Carlo 
simulations for all fillings, as long as U is less than the bandwidth and T is not much smaller than the crossover 
temperature Tx where renormalizcd-classical behavior sets in. Both quantum-critical and renormalized-classical 
behavior can occur in certain parameter ranges but the critical behavior of our approach is that of the O {n) model 
with n — > 00. p3| 

The main predictions of physical significance are as follows: 

1. The theory predicts a magnetic phase diagram where magnetic order persists away from half-filling but with 
completely suppressed ferromagnetism. 

2. In the renormalized classical regime above the zero-temperature phase transition, precursors of antiferromagnetic 
bands (shadow bands) appear in A(ki?,aj). These precursors occur when ^ > S,th (or losf < T). Between these 
precursors of antiferromagnetic bands a pseudogap appears at half-filling, so that the Fermi liquid quasiparticles 
are completely destroyed in a wide temperature range above the zero-temperature phase transition < T < Tx ■ 
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The upper critical dimension for this phenomenon is three. We stress the quahtative difference between the 
Hubbard bands and the precursors of antiferromagnetic bands and we predict that in two dimensions one may 
see both sets of bands simultaneously in certain parameter ranges. This prediction is consistent with the results 
of numerical simulations. [Q ||6^. We know of only one other analytic approach ^^ which leads to similar four 
peak structure in the spectral function. 

The zero temperature magnetic phase diagram is partly an open question because, despite the qualitative agreement 
with other analytical approaches, there is still an apparent contradiction with Monte Carlo simulations |64|. Our 
prediction of precursors of antiferromagnetic bands on the other hand is in agreement with Monte Carlo simulations. 
Neither this effect nor upper and lower Hubbard bands are observed in self-consistent schemes such as FLEX. This is 
because of inconsistent treatment of the vertex and self-energy corrections in this approximation, as we have explained 



in Sec.( Vl). However, if there was a Migdal theorem for spin fluctuations, it would be justifiable to neglect the vertex 
corrections and keep only the self-energy effects as is done in the FLEX approximation. The presence of precursors 
of antiferromagnetic bands in two-dimensions is then a clear case of qualitatively new Physics that would not appear 
if there was a Migdal theorem for spin fluctuations. The same is true for the Hubbard bands for large U > W in any 
dimension. 

We would like to state again clearly the nature of our critique of approximation schemes which are based on using 
Migdal's theorem for systems with electron-electron interactions. We do not imply that one does not need at all to 
take into account the feedback of the single-particle spectra on collective modes. The only point that we want to 
make here is that, based on sum rules and comparison with Monte Carlo data, we see that frequency and momentum 
dependent corrections to the self-energy and to the vertex often tend to cancel one another and that ignoring this 
leads to qualitatively incorrect results, in particular, with regards to the pseudogap. In this paper we were able to 
look only at the beginning of the renormalized classical regime when the pseudogap starts to form. The truly self- 
consistent treatment of the one-particle and two-particle properties in the pseudogap regime remains an open and very 
challenging problem. We hope that by extending our approach to the ordered state and looking at how the pseudogap 
starts to disappear as the temperature is raised, one can better understand how to develop a more self-consistent 
theory in the pseudogap regime. We now point out how our approach can be extended in other directions. 

As we mentioned in Sec.(M), the pseudogap and precursors of antiferromagnetic bands in the two-dimensional 
repulsive Hubbard model have interesting analogs in the attractive Hubbard model. In that model, one expects 
a pairing pseudogap and precursors of superconducting quasiparticle bands above Tc- At half- filling the negative 
and positive Hubbard models are mapped onto one another by a canonical transformation and the present theory is 
directly applicable to the attractive case. However, away from half-filling the mapping between the two models is more 
complicated and the microscopic theory requires additional sum-rule for pairing susceptibilities to find self-consistently 
the effective pairing interaction. This work is now in progress. 

The present approach can be also extended to stronger coupling U > W. Again the key idea would be to parame- 
terize the irreducible vertices, which have now to be frequency dependent, and then use the most important sum- rules 
to find the parameterization coefficients. This will, of course, require solving much more complicated self-consistent 
equations than in the present approach, but we believe that the problem still can be made tractable. 

Finally, we would like to make two comments about the magnetic and the pairing pseudogap in the context of 
high- Tc superconductors, based on the results of our studies. First, as was stressed in Ref. |Q, to understand clearly 
the physics of the single-particle pseudogap phenomena it is important to distinguish static short-range order from 
dynamical short-range order. The former is defined by a nearly Lorentzian form of the corresponding static structure 
factor S'(q) ex l/((q — Q)^ -I- C~^) (Q = (7r,7r) in magnetic case, Q = in the case of pairing), while the latter 
means only that the corresponding susceptibility x(qj 0) has such a Lorentzian form. A condition for the existence of 
the single particle pseudogap in the vicinity of a given phase transition is that the corresponding short-range order 
is quasi-static {i.e. ujsf ^ T) [p6| . Experimentally, one can measure directly the dynamical spin structure factor 
5(q, w), and then obtain the static structure factor through the integral S{q) = J S{q,uj)duj/{2TT). Even if the zero- 
frequency dynamical structure factor S'sp(q, 0) is very strongly peaked at q ^ Q it is possible that the static structure 
factor 5sp(q) is only weakly momentum dependent ||36|] . Thus in order to know whether one should expect to see the 
precursors of the antiferromagnetic bands and the corresponding pseudogap at a given doping and temperature it is 
necessary to obtain the static spin structure factor from the experimentally determined dynamical structure factor 
and then analyze its momentum dependence to see both if it is peaked and if it is quasi- two-dimensional. 

The second comment that we would like to make is that both the pairing and the magnetic single-particle pseudogap 
discussed above are an effect of low dimensionality and hence they exist as long as there is a large two-dimensional 
fluctuating regime before the real three-dimensional phase transition. In this context, a pairing pseudogap could exist 
on either side of optimal doping. |79| ] The much larger temperature range over which a pseudogap appears in the 
underdoped compounds suggests that, in addition to pairing fluctuations, other thermal fluctuations (charge, spin...) 
prohibit finite-temperature ordering. [ pO| An example of this occurs in the attractive Hubbard model where charge 
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fluctuations push the Kosterhtz-Thouless temperature to zero at half-fiUing, precisely where the crossover temperature 
to the pseudogap regime is largest. 
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APPENDIX A: SUM RULES, WARD IDENTITIES AND CONSISTENCY REQUIREMENTS 

In this Appendix, we recall well known constraints on many-body theory that follow from sum-rules and conservation 
laws and comment, wherever possible, on their physical meaning and on where commonly used approaches fail to 
satisfy these constraints. Although we come back on a detailed discussion of various theories in a later appendix, 
we find it useful to include some of this discussion here to motivate our approach. We consider in turn various 
results that would be satisfied by any exact solution of the many-body problem. They are all consequences of either 
anticommutation relations alone (Pauli principle) or of anticommutation relations and the Heisenberg equations of 
motion. We describe in turn: 1. The relation between self-energy and two-body correlation functions that embodies 
the details of the Hamiltonian; 2. Sum rules for one-particle properties; 3. Sum rules and constraints on two-particle 
properties, in particular /-sum rule and Ward identities that express conservation laws; 4. A few relations that are 
crucial in Fermi liquid theory, namely Luttinger's theorem and the forward scattering sum rule. 

1. Equations of motion and the relation between the self-energy S and two-particle properties. 

The self-energy (we always mean one-particle irreducible self-energy) is related to the potential energy, and hence 
to two-particle correlations through the expression Eq. (^) , which in the Kadanoff and Baym notation can be written 
as 

S,(l,T)a(T,l+)=C/(ntni). (Al) 

Here, the index with an overbar, 1, means that there is a sum over corresponding lattice positions and an integral 
over imaginary time. The notation l"*" means that the imaginary time implicit in 1 is ri -I- r/ where 77 is a positive 



infinitesimal number. Eq.(Al) is an important consistency requirement between self-energy and double occupancy in 
the Hubbard model that can easily be proven as follows. From the equations of motion for the single-particle Green's 
function Eq.(0) one finds, 



d 



Ga {vi,-V.j,Ti ~Tj) 



= <5.,,<5 (r, -T,)-U iTr [ct, {n) c.-. {n) C. {n) c], (r,)) ) . (A2) 

Using the short-hand notation in Eqs.(||) and (|j) and the definition of self-energy (Dyson's equation) the above 
equation is also written in the form. 

Go 1 (1,T) G„ (T,2) - <5(1 - 2) + E, (1,T) G, (T,2) . (A3) 

Comparing the last two equations, the well known relation Eq.([Al|) (or Eq.([44|)) between self-energy. Green's function 
and potential energy follows. 

So-called conserving P6[ approaches to the many-body problem violate the above consistency requirement Eq. (|44|) 
in the following sense. The right-hand side can be computed from the collective modes using the fluctuation-dissipation 
theorem. In conserving approximations, this gives a result that is different from what is computed directly from the 
left-hand side of the equation, namely from the self-energy and from the Green's function. In fact, all many-body 
approaches satisfy the above consistency requirement at best in an approximate way. However, it is a very important 
requirement and Eq.(H4) plays a key role in our discussion. Seen in Matsubara frequency, it is a sum rule, or an 
integral constraint that involves all frequencies, large and small. 
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2. Constraints on single- particle properties: 

The spectral weight A^ (k, uj) can be interpreted as a probabiHty of having an electron in a state (a, k, uj) and it 
satisfies the normalization sum rule 

1^ ^A, (k,c.) = ({ck.,cL}) = 1. (A4) 

Formally this is a consequence of the jump in the Green's function at t = 0, as can be seen from calculating 

G. (k,0-) - G. (k,0+) = 1 =Tj2{e^'^^' - e-'^-O H |^^i^l^ ^ /°° ^A^ (k,..) . (A5) 

To do perturbation theory directly for the Green's function to any finite order would require that the interaction U 
be small not only in comparison with the bandwidth W but also in comparison with the smallest Matsubara frequency 
iki = 27rr. Also, the direct perturbation series for the Green's function gives, after analytical continuation, poles 
of arbitrary high order located at the unperturbed energies. These high-order poles are inconsistent with the simple 
pole (or branch cut) structure of the Green's function predicted by the spectral representation. Furthermore, the 
high-order poles lead to a spectral weight that can be negative. |8l) The common way to get around these difficulties 
is to make approximations for the self-energy S instead and then calculate the Green's function using Dyson's equation 
Eq.®. 



It is interesting to note that to satisfy the constraint Eq.(A4), it suffices that I](k, i/c„), defined by Eq.rtq), has a 
finite limit as ifc„ -^ oo. More constraints on approximations for the self-energy may be found by continuing this 
line of thought. A systematic way of doing this is to do a high-frequency expansion for both the Matsubara Green's 



function and the self-energy and to find coefficients using sum-rules. The sum-rules that we need then are [p2| 

^^LoA, (k,c^) = ({[ck., {H - iiN)] , 4^}) = 6k - M + Un-. (A6) 



duj 



/CXD 7 

—oj-'A^ (k, oj) = (ek - A^)' + 2t/(ek - A*)"-. + U^n_, (A7) 

-oo 27r 

where Ua = n/2 since we are in the paramagnetic state. 

Using the spectral representation Eq.® one can easily see that the above sum rules give the coefficients of the 
high-frequency expansion of the Matsubara Green's function 



lim G^ (k,ifc„) = --+(-— ) / -^^A^ (k, w) -I- ( -- ) \ ^u?A„{V,uj) 



(A8) 



The self-energy has the same analytic properties as the Green's function. Using its high frequency expansion in the 
expression for the Green's function Eq.(P), one finds that the first term in Eq.([A^), leads to the requirement that the 
self-energy has a finite limit at i/c„ — ^ co. The second term fixes the value of this constant to the Hartree-Fock result, 
and the last and second-term combine to give the leading term in l/ik^ of the self-energy high-frequency expansion. 
In short, we find the result quoted in Eq.(|68|), namely 

lim S^ (k,zfc„) ^ Un_„ + ^'"'" ^] ~ ""'"^ + . . . (A9) 



The Kramers-Kronig relation for the self-energy 



Re [E- (k,.) - E« (k,oo)] . V f ^M^^ill:^ 



(AlO) 



and the high-frequency result Eq.( A9) imply the following sum-rule for the imaginary part of the self-energy 

- /" ^ Im [S]« ik,uj')] - U^n-, (1 - n_,) (All) 
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Important consequences of this equation are that for a given U the integrated imaginary part of the self-energy is 
independent of temperature and is increasing towards half-fiUing. The right-hand side of this equation is also a me asure 
the width of the single-particle excitation spectrum, as can be seen from the spectral weight moments Eqs.(A6) and 

(El), 



fj.'., (.,.,-[/ 



—ujA„(k,uj} 



U^n_, (1 



.) 



(A12) 



An important physical point is that the asymptotic behavior Eq^AE) is a necessary condition for the existence 
of upper and lower Hubbard bands, as has been explained in Sec. VI A. However, it is important to realize that 
it is not a sufficient condition. Indeed, the following paradox has been noticed in explicit calculations in infinite 
dimensions 174] , |ll[ . While ordinary second-order perturbation theory with bare Green functions Go reproduces 
correctly the appearance of the Hubbard bands with increasing [/, the perturbation theory with dressed Green 
function G = [Gq^ — S]^^ does not. The reason for this is that although the second-order expression for S^ (k,«A:„) 
in terms of full G does satisfy the asymptotics Eq.( A9), it sets in too late, namely for fc„ ^ U, instead of kn ^ W. 
The fact that the asymptotics shou ld start at kn ^ W even when U > W is a. non-trivial consequence of the Pauli 
principle, as explained in Sec. (VI A). Thus there are no Hubbard bands in any theory that uses self-consistent Green 
functions but neglects the frequency dependence of the vertex. This is an explicit example that illustrates what seems 
to be a more general phenomena: a calculation with dressed Green's functions but no frequency dependent vertex 
correction often gives worse results that the one done with bare Green's functions and a frequency independent vertex. 
We will see in the next subsection that this also happens in the calculation of the two-particle properties. Also, as we 
have argued in Sec. (VI), a similar situation occurs with the precursors of antifcrromagnetic bands in the renormalized 
classical regime in two-dimensions. 

Finally, we quote two more well known sum-rules that we will need. They involve the Fermi function / (uj) and the 
spectral weight. The first one follows from definition of G^ (k,r) and the spectral representation 



lim Go- (k,r) = 



J^ficu)A^ik,u;) 



'k(T 



Ckc 



n^o 



(A13) 



The quantity riko- is the distribution function. It is equal to the Fermi function only when the self-energy is frequency 
independent. The next result, that follows simply from the equations of motion. 



lim 



-E 



dG, (k,r) 
dr 



]V E / ^^-^ ('^) "^'^ ('^'^) = ]^ E ("^k - m) ni,a + U 



n^ni) 



(A14) 



k k " k 

is useful to show to what extent certain dressed-propagator approaches fail to satisfy the /—sum rule. 



3. Constraints on t^vo- particle properties 



For any one-band model, independently of the Hamiltonian, the Pauli principle (anticommutation relations) 

(n^,) = (n,;,) (A15) 

implies the following two simple identities: 

Uni^^±n^l) ) = n±2 (ni^riii) . (A16) 

The correlation functions on the left-hand side are equal-time and equal-position spin and charge correlation functions. 
The susceptibihties Xch (fz — Tj, t) , Xsp (fi — Tj, t) in Eqs.(p7|) and (16) are response functions for arbitrary (r^ — r^, t) 
so they must reduce to the above equal-time equal-position correlation functions when r^ = rvi and r = 0. This is one 
special case of the imaginary-time version of the fluctuation-dissipation theorem Eqs.(p^)(p7|). This translates into 
local-moment and local-charge sum-rules for the susceptibilities 



T 






J^Z^Z^Xsp (q, iqn) ^ 2 (n^n^) - 2 {m^ni) = n - 2 {m^ni) 



N 



(A17) 



T 






— 2^ 2^ Xch (q, iqn) = 2 (n^nij) + 2 (n^ni) -n =n + 2 {n^ni) - n 



N 



(A18) 



q iqn 
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where we have removed the i dependence of (ni^riii) using translational invariance. The right-hand side of the 

local-moment sum-rule is equal to ({S^) ) , while that of the local-charge sum rule is equal to (p^) — n^ 

If arbitrary sets of diagrams are summed, nothing can prevent the right-hand side from taking unphysical values. 
For example, the Pauli principle ma y be violat ed, i. e. {n-^n^) j^ (n|) . To see this, notice that when the Pauli principle 
is satisfied, our two sum rules Eqs.(|A17) and (A18) lead to 



T 



T7 X! X! ^^^P (^' *9") + ^"f' (^' ^^")] = 2n - n^ 



TV 



(A19) 



q iqn 



It is easy to check that well known approaches to the many-body problem, such as RPA, violate this basic requirement. 
Indeed, the ordinary RPA expressions for spin and charge are 



RPA 

A-sp 



ii) 



Xo 



rXo 



(A20) 
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l + fxo 



(A21) 



where 



T 



Xo(g)--2-5]GW(fc)G(°'(fc + g). 



(A22) 



That RPA does not satisfy the sum rule Eq.(Al9) already to second order in U can be easily seen by expanding the 
denominators. 

To satisfy the Mermin-Wagner theorem, approximate theories must also prevent {n^ni) from taking unphysical 
values. This quantity is positive and bounded by its value for U — oo and its value for non- interacting systems, 
namely < {n-^n^) < n^/4. Hence, the right-hand side of the local-moment sum-rule Eq.(Al7) is contained in the 



interval [n, i 



Any theory that prevents the right-hand side of the local-moment sum rule from taking infinite 



values satisfies the Mermin-Wagner theorem. 



Proof: Near a magnetic phase transition, the zero Matsubara-frequency component of the spin susceptibility takes 
the Ornstein-Zernicke form 



Xsp(q + Q,0) 



1 



-e- 



(A23) 



where q is measured with respect to the ordering wave vector Q and where ^^ is the square of the correlation 
length. Near its maximum, the above susceptibility is of order ^^ while all finite Matsubara-frequency compo- 
nents at the ordering wave vector are at most of order 1/ {2ttT) which is much smaller than ^^. Hence, one 
can keep only the zero-Matsubara frequency contribution on the left-hand side of the local-moment sum rule 
Eq.(|A17|) obtaining 



T 



r _d^ 

J (27t)' 



(27r)'' 9^ 



C 



(A24) 



where C contains non-zero Matsubara frequency contributions as well as n — 2 (nini) . Since C is finite, this 
means that in two dimensions (d = 2), it is impossible to have ^"^ = on the left-hand side otherwise the 
integral would diverge logarithmically. 



Finally, the /-sum rule on spin and charge susceptibilities follows as usual from the fact that the Hamiltonian 



conserves particle number. Computing 



Pq 



dp- 



and 



T = 



dS- 



It=0 



one obtains for either charge or spin 
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(A25) 
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As can be seen from the spectral representations of spin and charge susceptibihties, Eq.(gO|), the quantity that obeys 
the /—sum rule is the coefficient of the leading term in the l/g^ high-frequency expansion of the susceptibilities. 

The single-particle energies ek entering explicitly the right-hand side of the /—sum rule are independent of inter- 
actions, so interactions influence the /—sum rule only very weakly through the nto-. In fact, in a continuum et oa k^ 
so nkcr enters only in the form X^ku ""ko- — n . In this case, the right-hand side of the /—sum rule is proportional to 
q^n and hence is independent of interactions. On a lattice however, the energies cannot in general be taken out of the 
sum and interactions influence the value of the right-hand side, but only through the fact that nko- differs from the 
non-interacting Fermi function /ko-. At strong-coupling, where the self-energy is strongly frequency dependent, this 
difference between riko- and /ko- becomes important. But from weak to intermediate coupling, calculations where /ko- 
appears on the right-hand side should be good approximations. In the explicit examples that we have treated, the U 
dependence of the /—sum rule becomes important only close to half-filling and for f/ > 4, signaling the breakdown of 
approximations based on frequency-independent self-energies. 

While RPA-like theories that use /ko- instead of nko- violate only weakly the /— sum rule in the weak to interme- 
diate coupling regime, self-consistent theories that use frequency-dependent self-energies but no frequency-dependent 
vertices violate conservations laws in general, and the /— sum rule in particular, in a much more dramatic way. The 
point is that susceptibilities with a dressed bubble, xrpa = Xo/(l ~ o^^o)' ^'"^ ^^^ approximations because they 
have the following properties, for any value of U 



Xrpa{(1 = O.iQn 7^ 0) 7^ 



(A26) 
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(A27) 



The first of these equations explicitly violates the Ward identity, Eq.(A31) below, at all frequencies, including 
small non-zero ones, since at zer o wa ve vector we should have x( q = 0,iqn ^ 0) = for all frequencies except 
zero. The second equation, Eq.(A27) violates the /-sum rule Eq.(A2J;) at all wave vectors, by a constant term 
4[/ {{ni) (ni) — (7i|7i|)) which in practical calculations, say ai U — A, is of the same order as the first term, which is 
the only one that should be there according to the /-sum rule. 



Proof: Eqs.( |A26| ) and (A27) are proven as follows. Consider the standard RPA expression but with dressed bubbles 



XRPA 



Xo/(l- -^Xo)- 



(A28) 



Using the spectral representation for the Green's function and inversion symmetry in the Brillouin zone one 
finds. 



Xo 



(q, *9n) = I E / ^ / ^^ (1^'^) ^ ("^ + q^ 



,, (a;-c^O(/(c^O-/M) 

UJ ) 5 . 

{co-c'f+ql 



(A29) 



When the bubble is not dressed, the spectral weights are delta functions so that at q = the susceptibility 
would vanish for all non-zero values of (/„, as required by the Ward identity. However, here because the spectral 
weight has a width and because the integrand is even and positive, then the integral will not vanish, resulting 
in the first anomaly Eq.(A26) we mention. To prove the second equation, Eq.(A27), it suffices to reme mbe r 
from the spectral representation of the susceptibility Eq.(|2(]|) and the derivation of the /—sum rule Eqs.( [A25| ) 
that we are looking for the coefficient of the l/q^ term in the high-frequency expansion. Given the RPA form, 
Eq.(A2S), only the numerator contributes to this limit. One obtains, for the coefficient of the 1/g^ term, 



N 



E 



duj 

2^ 



du' 

"2^ 



A {k,oj) A{k + q,c^') {io - u;') (/ {u;') - / {u;)) 



(A30) 



from which Eq.(A27) follows using the sum rules for occupation number Eq.(A13) and for energy (A14) 



Conservation laws have general consequences not only on equal-time correlation functions, as in the /—sum rule 
above, but also on time-dependent correlation functions. For example, from the Hcisenberg equations of motion and 
anti-commutation relations, follow the Ward identities lEsI 



38 



Y.Y. E ( ^ + (^k+q - fk) ) (Trcl^ (r) a^Ck+q. (r) 4,+^^, (n) a'^Ck'.' (ra) 



k rT=±lcr' = ±l ^ 

= <5(r-Ti) ^ a'^G,,(k',r2-T)-^(T-r2) ^ a'^G,- (k' + q, r - n) (A31) 

ct' = ±1 (t'=±1 

where ^ = for charge, and i — 1 for spin. The /—sum rule above, Eq.( |A25| ), foUows from the above identity by 
simply taking ti = t^ , summing over k' and subtracting the two results for t ^ t^ and t —> t^ . 

We have seen in this section that there are strong cancelations for two-particle properties between the frequency 
dependence of self-energy and that of the vertex corrections, so that putting a frequency dependence in only one of 
them is a bad approximation. We have adopted the Kadanoff-Baym formalism in the main text since it can be used 
as a guide to make approximations that satisfy conservation laws. 

4. When there is a Fermi surface 

When perturbation theory converges (no phase transition) then at zero temperature T = the imaginary part of 
the self-energy vanishes, EJ^ (k,a; = 0) = 0, for all k values and the Fermi surface defined by 

ek - A* - S;, (k,c^ = 0) = (A32) 

encloses a volume that is equal to the volume enclosed by non-interacting particles 

1 ^ (^ - fk - S; (k,0)) = l-Y.0{f,,-e^)= n,. (A33) 

k k 

This is the content of Luttinger's theorem. |2^ |8^. It implies that there is a strong cancelation between the change 
of the chemical potential and the change of the self-energy on the Fermi surface. In particular, when EJ^ (kF,0) does 
not depend on k or on the direction of k^? (infinite D Hubbard model, electron gas) the change in (fi — (Iq) is exactly 
canceled by E^ (ki?,0) 



/i-/io = S;(kF,0). (A34) 



Luttinger's theorem is satisfied when 



lim f^^^J^^Ga{k,iiy)diydk = 0. (A35) 

Any theory that calculates its self-energy from a functional derivative of the Luttinger-Ward functional E = S^[G]/SG 
will satisfy Luttinger's theorem |2q ] |p3| . The latter procedure requires self-consistent determination of the self-energy 
as a function of momentum and frequency E^- (k,ifc„) and is usually quite computationally involved. However, even 
when this procedure to calculate the self-energy is not followed, it turns out to be rather easy to satisfy this theorem 
to an excellent degree of approximation in the weak to intermediate coupling regime. The reason for this is that 
any frequency-independent self-energy will preserve Luttinger's theorem and weak frequency dependence will not 
cause great harm. For the electron gas, Luttinger |28] suggests a way to build a perturbation theory in terms of 
non-interacting Green's functions which allows to satisfy Luttinger's theorem to very good accuracy. The trick is that 
the chemical potential for the interacting electrons /i should always enter the calculations in the form of the difference 
with the shift of the self-energy on the Fermi surface Go = l/[ikn — ek -I- (/^ — E^ (kF,0))]. The "non-interacting" 
Green's function Go in this formalism is the Green's function of some effective non-interacting system and, in general, 
it is different from both 1/ (ifc„ — Ck + m) ^^^id 1/ {i^n ^ Ck + Mo)- However, when T ^ Luttinger's theorem requires 
that (/i — EJj (ki?,0)) —f fiQ and one can approximate Go by the Green's function for a non-interacting system of the 
same density Go = 1/ {ikn — Ck + Mo)- In practice, one can also have a phase transition (or crossover) at a finite 
temperature T^ (Tx)- In these cases Luttinger's theorem is satisfied only approximately since the zero-temperature 
limit cannot be reached without a breakdown of perturbation theory. Then the relevant question is how well it is 



satisfied at Tc [Tx)- (See also Sec.([IIB 2) for a discussion of Luttinger's theorem in our approach) 

When Luttinger's theorem holds, one can usually develop a Landau Fermi liquid theory. In this approach, the Pauli 
principle is implemented only for momentum states near the Fermi surface by imposing the forward scattering sum 
rule. This sum rule, in two dimensions, reads 
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E 



f; 



Ff 



1 + f; 1 + Ff 



= 



(A36) 



where F^" and F° are the symmetric and antisymmetric Landau parameters expanded on the e~ basis instead of 
the Legendre polynomial basis. Recent renormalization group analysis has however claimed |84[| that the forward 
scattering sum rule comes from an inaccurate use of crossing symmetry and is not the proper way to enforce the Pauli 
principle. Most approaches to the many-body problem disregard this sum rule anyway, in the same way that they 
disregard the local Pauli principle. 



APPENDIX B: PROOFS OF VARIOUS FORMAL RESULTS 



In this appendix, we give the proofs of various relations mentioned in Sees. (Ill) and (IIIBS) 



1. The general expression for the self-energy Eq.(p7|) can be obtained as follows. Use the equations of motion and 
the definition of the self-energy Eqs.(A2)(|A3|) which in the present notation give 



S, (1, 1) G, (1, 2) = ~U (Tr [^+, (1++) v.-. (1+) v. (1) ^t (2)] ) 



(Bl) 



= -U 



|^-C_.(l.l^C,(l,2) 



(B2) 



Substituting the equation for the three-point susceptibility (collective modes) Eq.(|2q) in this last equation and 
multiplying on both sides by G~^ proves [E^ the expression Eq.(p7|) for the self-energy. 

2. We now show that our approach satisfies the consistency requirement between single-particle propertie s and 
collective modes in the form of Eq. ([4q) . Using our expression Eq.(|4q) for S^^^ and the definition of xo Eq.( [A22| ) 
we obtain 



1™ ^ E ^^'^ (^) G^"^ (^) ^''""^ = ^"- - T^ E [UspXspiq) + Ur^hXchiq)] ^ 



:*o- TV 



(B3) 



Using 



Xsp{q) - xq (q) = -|^xo (q) Xsp{q) 



(B4) 



xo{q) -Xch{q) 



Uch 



xoiq)xchiq) 



(B5) 



and the local moment Eq.(BS) and local charge Eq.(37) sum rules proves the result. The result is also obvious 



if we follow the steps in the first part of this Appendix to deduce the self-energy expression Eq.(31) using the 
collective mode equation Eq.(p(l) adapted to our approximation. 



APPENDIX C: ANSATZ FOR RELATION BETWEEN Usp AND {N^Ni). 



Using the present notation and formalism, we now give a p hysical derivation of Eq.(|o|) that is equivalent to the 
one already given using the equations of motion approach |29|| . (The latter derivation was inspired by the local field 
approximation of Singwi et al. pl| ) . Since our considerations on collective modes are independent of the precise value 
of the interaction U , we do have to use the equations of motion, or the equivalent, to feed that information back in 
the definition of irreducible vertices. The two irreducible vertices that we need are in principle calculable from 



'(5 (1-3) (5 (2 -4) (5 (2-1+) = 



6^„ (1, 2) _ 5 [S. (1, 1) G. (1, 2) G-i (2, 2)] 



<5G,,(3,4) 



^G,'(3,4) 



(CI) 
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The rewriting on the right-h and side has been done to take advantage of the fact that in the Hubbard model, the 
equations of motion (see Eqs.(A2)(A3)) give us the product So- (l, l) Go- (l, 2) as the highly local four field correlation 

function — C/^7V [^-o (l^^)'0-cr (l^)'0cr {^)4't (2)]) ■ Ordinary RPA amounts to a Hartree-Fock factoring of this 
correlation function. Pursuing the philosophy that the minimum number of approximations should be done on local 
correlation functions, we do this factoring in such a way that it becomes exact when all points are identical, namely 
when 2 = 1^. In other words, we write 

- U (T. [V.+0 (1^+) ^-. (1+) ^. (1) i^t m)- ^ (l7(l)|"^/[V - (1' 1^) G^ (1' 2) (C2) 

All quantities are evaluated as functionals of G up to this point. Wc can now evaluate the functional derivative 



<5Go' (3,4) 



U 



("t(i)"-i(i)) 

(nT(l))(«i(l)) 



G_o(l,l+)'5(l-2) 



<5Go'(3,4) 



(C3) 



U 



("T(l)"i(l)) 
!"T(l)><»i(l)> 



5Go' (3,4) 



G_o(l,l+)'5(l 



(nT(l)ni(l)) <5G_o(l,l+) 
(n^ (1)) (n, (1)) 5Go'(3,4) ^ ^ 



(C4) 



The functional derivatives are now evaluated for the actual equilibrium value of G. Hence, we can use rotational 
invariance to argue that the first term is independent of a and tr' whereas the last one is proportional to 5-cr.a' ■ Since 



that 



Fii, only this last term proportional to S 
'5G_o(l,l+) 



/ contributes to U.... To obtain this term, it suffices to note 



-^ sp 



<^Go' (3,4) 



= (5_o,a"^(l-3)(5(4-l+) 



(C5) 



and we obtain the desired result Eq.(40) for U, 



^ sp- 



APPENDIX D: REAL-FREQUENCY ANALYSIS OF THE SELF-ENERGY AND FERMI LIQUID LIMIT 

It is instructive to recover the two-dimensional result for precursors of antiferromagnetic bands using the real- 
frequency formalism since it also clarifies the limit in which the Fermi liquid result is recovered. Again we neglect 
the contribution of charge fluctuations. Starting from our expression for the self-energy Eq.(^), one uses the spectral 
representation for the susceptibility and for G'-"-' . The Matsubara frequency sums can be then done and the result is 
trivially continued to real frequencies, pq] One obtains, for the contribution of classical and quantum spin fluctuations 
to the self-energy in d dimensions 



£^(k,u;) 



UU, 



df^q 

(2^ 



dJ 



n(^') + /(ek+q)] 



x'L (q- ^') 



< + 11-I + UJ' - (ek+q - A*o) 



(Dl) 



where /iq = at half-filling in the nearest-neighbor model and where / is, as usual, the Fermi function, while 
n {uj) = (e^'^ — l) is the Bose-Einstein distribution. To analyze this result in various limiting cases we need to 
know more about the frequency dependence of the spin susceptibility. When the antiferromagnetic correlation length 
is large, the zero- frequency result Eq.(pO) mentioned above can be generalized to 






1 



1 4- q2^2 _ iu}/uJsF 



where, ujsf = D /^^ is the characteristic spin relaxation frequency. In the notation of Ref. 
diffusion constant D is defined by 



1 = !° 



(D2) 
the microscopic 

(D3) 



with the microscopic relaxation time, 
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To 



dxo (Qd,w) 



xo [Qd) 



diuj 



(D4) 



j=0 



This relaxation-time is non-zero in models where the Fermi surface intersects the magnetic Brillouin zone. Clearly, 
the frequency dependence of x^(q+ Qd;^) is on a scale ujsf = D/^"^. The l/w decrease of x'sp at high-frequency 
is not enough to ensure that the real frequency version of the local-moment sum rule is satisfied and the simplest 
way to cure this problem is to introduce Kq] a high-frequency cutoff ^cut- The large correlation length makes the 
characteristic energy of the spin fluctuations ujsf a small number (critical slowing down). We consider in turn two 
limiting cases. [Q The Fermi-liquid regime appears for lose ^ T and the non-Fermi liquid regime in the opposite 
(renormalized classical) regime ojsf ^ T. 



1. Fermi liquid and nested Fermi liquid regime ujsf ^ T 

Perhaps the best known characteristic of a Fermi liquid is that S"^ (kj?, lj; T = 0) ex lj^ and E"^ (k^, lj = 0; T) cx 
T^ . To recover this result in the regime ujsf 3> T far from phase transitions, we start from the above expression 
Eq.(Dl) for the self-energy to obtain 



S"^(kF,^) 



UU..„ 1 



Jd-l 



9± 



duj 



[n {uj') + f{u + uj')] x% {q^,q\\ {q^,^F,uj, uj') ; J) 

TT /^ \ II / 



(D5) 



4 2vF J (27r)" 
where g|| , the component of q parallel to the Fermi momentum hp , is obtained from the solution of the equation 

£k+q - fJ.O = UJ + Uj' (D6) 

The key to understanding the Fermi liquid vs non-Fermi liquid regime is in the relative width in frequency of 
x'sp (Qi^O /^' I's the width of the combined Bose and Fermi functions. In general, the function n (w') + f {uj + uj') de- 
pends on uj' on a scale Max {uj, T) while far from a phase transition, the explicit frequency dependence of x'' (q,i^') /uj' 
is on a scale ujsf ^ Ef ^ T . Hence, in this case we can assume that x'sp (l;'^') 1'^' is a constant in the frequency 
range over which n {uj') + f {uj + uj') differs from zero. Also, since x'sp {^^^') /w' depends on wave vector q over a scale 
of order qF, one can neglect the uj + uj' dependence of g|| obtained from Eq.(D6). Hence, we can approximate our 
expression Eq. (p5|) for E"^ by 



S"^(kf,c.)~- 
where the substitution x 



UUspA{]<iF) f duj' 



[n {uj') + f{uj + uj')] uj' 



UUspA{kF) 
4 2vf J tt ''''''' '" " 4 Avf 

e^'^ allowed the integral to be done exactly and where 



{nTY 



A{kF) 



f d'^'^q^ Xl'p(g±,'7||('?±,kf,0,0);cj') 



{2n) 



d-l 



lim 



(D7) 



(D8) 



In general, A depends on the orientation of the Fe rmi wave vector, kj?, because it determines the choice of parallel 
and perpendicular axis q±,q\\. The above result Eq.(D7) for E"^ is the well known Fermi liquid result. 

There are known corrections to the Fermi liquid self-energy that come from the non-analytic uj' jvpq behavior of 
x" (q,a;') juj' near the ferromagnetic (zone center) wave vector. In three dimensions ^8| this non-analyticity leads to 
subdominant w^lnw corrections, while in two dimensions it leads to the dominant oj^lnw behavior. |^ |Q In the 
case under consideration, the antiferromagnetic contribution has a larger prefactor. Even when it dominates however, 
it can also lead to non-analyticities in the case of a nested Fermi surface. Indeed, we note that 



Imx^(Qd,c^) - ^iVd(|)tanh (^ 



(D9) 



In two dimensions, the logarithmic divergence of the densit y of states Nd{%) at the van Hove singularity makes the 
zero- frequency limit of the microscopic relaxation time Eq.(D4) ill-defined, because of the logarithmic divergence at 
UJ = 0. However, this leads only to logarithmic corrections. If we drop logarithmic dependencies then for uj < T 
one has dxa {Qd,uj) ldiuj\^^rj. ^ 1/T and this l/T dependence of dxo {Q,diUj) / diuj\^^^ changes the temperature 
dependence of Y."^ (ki?, 0) from T^ to T as discussed in the "Nested Fermi Liquid" approach. [pi| 
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2. Non- Fermi liquid regime lusf "C T 



Near an antiferromagnetic phase transition, the spin-fluctuation energy becomes much smaher than temperature. 
This is the renormahzed classical regime. The condition ujsf ^ T means that Xsp (^-L, 9||;^') is peaked over a 
frequency interval w' <C T much narrower than the interval w' ~ T over which n (tu') + f [lu + lu') changes. This 
situation is the opposite of that encountered in the Fermi liquid regime. To evaluate E"-^ Eq.(pl)5|) the Fermi factor 
can now be neglected compared with the classical limit of the Bose factor, T/cu'. Then the dominant contribution to 
T,"^ {kp,uj) is from classical spin fluctuations T J ^-^x'sp = ^Xsp — Sgp as we see below. More specifically, we take 
into account that the integral is peaked near Q = (tt, tt) and measure wave vector with respect to the zone center. For 
simplicity we consi der b elow the half-filled case /Lto = 0. Then, with t he h elp of ek+q+Q = ^£k+q we approximate the 
equation for qn Eq.(D6) by vpqu = — {lo + u>') . This gives us for Eq.(D5) the approximation 



S"«(k^,c.)«- 



Jd-l, 



UUsp 1 
4 2vF J (27r) 



sd-l 



duj' T ,, 

— -,Xsp\<l^:q\\ 



Vp 



-;^ 



(DIO) 



The dependence of x'sp oii ^' through ^n = — (a; + uj') /vp may be neglected because ^n appears only in the 
combination i^~^ ~^ 1± ~^ lf\) ^^'^ ^^ the regime usf <^ T we have to'/vp < ujsf/vp ^ D^~'^ /vp <C C^^. The latter 
inequality is generically satisfied when S^~^ <C 1. Using 
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duo' 1 



iXsp 9-L>9|| 



Vp 



Tx'sp 1 <i-i-.q\\ = ;«9ii 

^ \ Vp 
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\vf) 



(Dll) 
(D12) 



the above equation Eq. (|D10|) for Y."^ {kp 

time formalism Eq.(|5g). As we saw in Sec.(VAl 

dominant and leads to limr^o E"^ (kp, 0) 



reduces precisely to the classical contribution found using imaginary- 
when the condition ^ > ^^^j is satisfied, then this contribution is 



APPENDIX E: EXPANDED DISCUSSION OF OTHER APPROACHES 



This appendix expands on Sec. (VIII) to discuss in detail various theories, explaining the advantages and disadvan- 
tages of each in the context of the sets of constraints described in Appendices (A 2) and ([A3|). 



1. Paramagnon theories 



In standard paramagnon theories, [ 
bare or dressed Green's functions. T 
variety of paramagnon theories. They are largely 
We concentrate our discussion on recent versions 



p2] the spin and charge fluctuations are computed by RPA, using either 

ren the fluctuations are fed back in the self-energy. In fact there is a whole 

henomenological. The reader is referred to Ref. Eq] for a review. 

i2] of the so-called Berk-Schrieffer formula P3| . In this approach, 



infinite subsets of diagrams are summed and bare propagators are used in the calculation of both the susceptibilities 
and the self-energy, the latter being given by 



^BS 



(k) = Un^, + J^J2 [{^Uxff^iq) 2Uxoiq)) 



-Uxfnq)]G'Ak + q). 



(El) 



The RPA spin and charge susceptibilities have been defined in Eqs.( A20| )(A21). Comparing with our self-energy 
formula Eq.(p^), it is clear that here there is no vertex correction. In addition, the factor of three in front of the spin 
susceptibility in Eq.(El) is supposed to take into account the presence of both longitudinal and transverse spin waves 
and the subtracted term is to avoid double-counting the term of order U"^. 

We can now see the advantages and disadvantages of this approach. First, note that the susceptibilities entering the 
Berk-Schrieffer formula are the RPA ones. As we saw in Appendix(H), these fail to satisfy both the local Pauli principle 
and the Mermin- Wagner theorem. Hence, spurious phase transitions will influence the self-energy in uncontrollable 
ways. The collective modes do however satisfy conservation laws since they are obtained with bare vertices and Green's 
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functions containing a constant self-energy. The /-sum rule Eqs.(A25) then is satisfied without renormalization of the 
distribution function nk because the zeroth order self-energy is constant. This is all in agreement with the definition 
of a conserving approximation for the collective modes. 

The high-energy asymptotics of the self-energy sets in at the correct energy scale kn > W in this approach, but the 
second term of the large-frequency asymptotics is incorrect. Indeed, at large values of zA:„, 



lim E^* (fc) == C/n_„ + 

ikn — >-CXj 



Y^ 



U T 

Aikn N ^-^ 
1 



[sc/xf/-"!?) + t/x^r (<z) - 2L/xo(g)] + 



and the sums can be evaluated as follows using the fluctuation-dissipation theorem 



(E2) 



(E3) 
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(q) = 2 (nt?it) + 2 (n|n|) - n^ 



(E4) 
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(E5) 



The correlators on the right-hand side take their RPA value so they do not satisfy the Pauli principle, i.e. {n-\ni) ^ 
(n|) . Taking these results together we have 



2 r 



lim Sf ^ (fc) = L/n_, + -^ 

ikn-^oc IrCn 



2("T"T> - ("T'^i) 



(E6) 



This does not gives the correct asymptotic behavior Eq.fcS) even if the Pauli principle {n^n^) = {n^) were satisfied, 
because {n-\n^) depends on the inte rac tion U . 

The paramagnon self-energy Eq.(El) also does not satisfy the consistency requirement Eq.(p^ between self-energy 
and collective modes imposed by the equations of motion. To see this we first note that 

1™ Ij E ^^' (^) G^"^ (^) ^"^'"^ = f^-- - ?^ E [3X^^(9) + Xr^(<z) - 2xo(g)] Xo (g) (E7) 



.0- N 



8 TV 



Using this expression in the sum-rule Eq. (Bq) which relates one and two-particle correlators and expanding both sides 
of this sum- rule in powers of C/, one finds that it is satisfied only up to order U^ . On the other hand, if one replaces 
2>Xsp — 2xo in Eq.(El) by Xsp, the sum-rule Eq.(Bq) is satisfied to all orders in U. In our opinion, the problem of 
enforcing rotational invariance in approximate theories is highly non-trivial and cannot be solved simply by adding 
factor of 3 in front of Xsp and then subtr acting 2 yo to avoid double counting. For more detailed discussions see Ref. 
^0| and the comments at the end of Sec.( [IIB 2 ). 

Luttinger's theorem is trivially satisfied if the occupation number is calculated with the initial constant self-energy 
since it gets absorbed in the chemical potential. If the occupation number is calculated with the Green's function 
that contains the Berk-Schrieffer self-energy then Luttinger's theorem is in general violated. It is advisable to use a 
new chemical potential. 



2. Conserving approximations (FLEX) 



In the conserving approximation schemes ||2^], one takes any physically motivated subset of skeleton diagrams 
to define a Luttinger-Ward functional $. Skeleton diagrams contain fully dressed Green's functions and no self- 
energy insertions. This functional is functionally differentiated to generate a self-energy that is then calculated 
self-consistently since it appears implicitly in the Green's functions used in the original set of diagrams. A further 
functional differentiation allows one to calculate the irreducible vertices necessary to obtain the collective modes in a 
way that preserves Ward identities. If one uses for the free energy the formula 



\nZ = Tr [In (-G)] + Tr (SG) - $ 



(E8) 



44 



then one obtains thermodynamic consistency in the sense that thermodynamic quantities obtained by derivatives of 
the free energy are identical to quantities computed directly from the single-particle Green's function. For example, 
particle number can be obtained either from a trace of the Green's function or from a chemical potential derivative of 
the free energy. In this scheme, Luttinger's theorem is satisfied as long as perturbation theory converges since then 
any initial guess for the Luttinger-Ward functional will satisfy Luttinger's theorem. 

FLEX refers to a particular choice of diagrams for $. This choice leads to the following self-consistent expression 
for the self-energy 

Ef ^ (fc) = [/n_. + 7§ E [(.^Uxff^iq) - 2UUq)) + Ux^f^^iq)] G^k + q). (E9) 

This expression for the self-energy does not contain vertex corrections, despite the fact that, contrary to the electron- 



phonon case, Migdal's theorem does not apply here. We have explained in detail in Sec.(VIB) why this may lead 
to qualitatively wrong results, such as the absence of precursors of antiferromagnetic bands and of the pseudogap in 
A(kp,Lu) in two dimensions. 

Another drawback of this approach is that it does not satisfy the Pauli principle in any form, either local or through 
crossing symmetry [^. Indeed, one would need to include all exchange diagrams to satisfy it. In practice this is 
never done. In the same way that there is nothing to constrain the value of (n-^n-^) obtained by the fluctuation- 
dissipation theorem to be equal to (fi-i), there is nothing to explicitly constrain the value of (n|n|). Nevertheless, the 
Mermin- Wagner theorem is believed to be satisfied in FLEX because the feedback through the self-energy tends to 
prevent the divergence of fluctuations in low dimension. [p4|, l38|. Physically however, this seems to be an artificial 
way of satisfying the Mermin- Wagner theorem since this theorem should be valid even in localized spin systems where 
single-particle properties are negligibly influenced by thermal fluctuations. We also point out that the proof of the 
Mermin- Wagner theorem in n — > oo models implies that the finite temperature phase transition in two dimensions 
is not simply removed by thermal fluctuations, but that it is replaced by a crossover to the renormalized classical 
regime with exponentially growing susceptibility. The fact that the conserving susceptibility in FLEX does not show 
such behavior |3^ means that FLEX is actually inconsistent with the generic phase space arguments responsible 
for the absence of finite-temperature phase transition in two dimensions. The case of one dimension also suggests 
that collective modes by themselves should suffice to guarantee the Mermin- Wagner theorem without feedback on 
single-particle properties. Indeed, in one dimension one shows by diagrammatic methods (parquet summation or 
renormalization) that the zero-temperature phase transition is prohibited at the two-particle level even without self- 
energy effects [g . 

Although, the second-order diagram is included correctly in FLEX, it does not have the correct coefficient in the 
1/ikn expansion of the self-energy. More i mportantly, the high-frequency behavior sets-in too late to give the Hubbard 



bands, as we have explained in Sec.( VIB ). We have also seen a case where FLEX, as judged from comparisons with 
Monte Carlo simulations (Fig. (la) of Ref. pQ]), does not reproduce the results of second-order perturbation theory 
even when it is a good low-energy approximation. 

One of the inconsistencies of conserving approximations that is seldom realized, is that the self-energy is inconsistent 
with the collective modes. In other words, the consistency formula Eq.(p4|) is not satisfied in the following sense. 
The explicit calculation of SG leads to an estimate of U (n|n|) that differs from the one obtained by applying the 
fluctuation-dissipation theorem to the conserving spin and charge susceptibilities. 

3. Pseudo-potential parquet approach 

In the parquet approach, one enforces complete antisymmetry of the four point function by writing down fully 
crossing-symmetric equations for these. There are three irreducible vertices, namely one for the particle-particle 
channel, and one for each of the two particle-hole channels. They obey the so-called parquet equations. |£^ The Green's 
functions are dressed by a self-energy which itself contains the four point function. In this way, self-consistency between 
one-particle and two-particle quantities is built-in. Solutions are possible for the one-impurity problem |96|] and in 
one-dimension Iq]. However, to solve the parquet equations in higher dimension with presently available computing 
power is impossible. Bickers et al. ||25|| ||53| have formulated the parquet equations as a systematic improvement over 
FLEX and have devised a way to do practical calculations by introducing so-called pseudo-potentials. Since the main 
computational difficulty is in keeping the full momentum and frequency dependence of the four point functions entering 
the calcu lation of the self-energy, this is where the various fluctuations channels are approximated by RPA-like forms 



Eq.(A2S) but with fully dressed propagators and an effective interaction (pseudo-potential) instead of U . A different 
strategy is under development. P4] The criticism of the present section applies only to the current pseudo-potential 
parquet approach. [2^ |5^ 
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It can be seen that one drawback of this approach at the physical level is that the use of constant effective 
interactions with dressed single-particle propagators means that the fluct uatio ns used in the calculation of the self- 
energy do not satisfy conservation laws, as we just demonstrated in Sec(A2). Furthermore, the pseudopotentials 
are determined by asking that the susceptibilities extracted from the four-point functions in the parquet equations 
match the corresponding RPA-pseudo-potcntial susceptibility at only one wave vector and frequency. The choice of 
this matching point is arbitrary: Should the match be done for the typical, the average, or the maximal value of the 
susceptibility in the Brillouin zone? 

As we have seen in Sec. (VI), even if the expression for the self-energy in this approach explicitly has the second- 
order perturbation theory diagram in it, this is not sufficient to ensure that the correct high frequency asymptotic 
behavior starts at the appropriate frequency scale ikn ~ W. Nevertheless, in many cases the results of the calculations 
performed with this approach are not so different from second-order perturbation theory, as can be seen from Fig.(l) 
ofRcf.( |§). 

Going rapidly through the rest of our list of properties, we see that the consistency requirement 
Scr (l,T) Ga (l, 1^) = U {n^ni) is at least approximately built-in by construction. Concerning the local-moment sum- 
rule and the Mermin- Wagner theorem, it has been shown that the so-called "basic" parquet equations should have the 
same critical behavior as the leading term in the 1/A^ expansion, [ p7[ and hence should satisfy the Mermin- Wagner 
theorem. [Q The pseudo-potentials should not affect the self-consistency necessary to satisfy the Mermin- Wagner 
theorem but the fact that they are matched at a single point might introduce difficulties, especially if the wave- vector 
at which Xsp becomes unstable is unknown from the start. As far as the Pauli principle is concerned, it should be 
at least approximately satisfied both locally and in momentum space. Nothing however in the approach enforces 
conservation laws. 



4. Present approach 



The role of the above sum-rules in our approach has been discussed in detail in the main text. Here we will discuss 
only a few additional points. 

If we concentrate on the q = properties, our spin and charge correlations behave as a special case of the "local 
Fermi liquid" defined in Ref. |g^. A "local Fermi liquid" is a description of q = properties that applies when 
the self-energy, and consequently irreducible vertices, depend only on frequency, not on momentum. In a local Fermi 
liquid there are only two Landau parameters, which in our case are Fq = —UgpXo (0^, 0) /2 and Fq — UchXo (0^, 0) /2. 
Unitarity and the forward scattering sum rule, if valid, imply that there is no ferromagnetism in the repulsive case 
|]98| , as we have found. One can check explicitly that the forward scattering sum rule is satisfied to within about 
15% in our usual Monte Carlo parameter range. However, as discussed in Appendix (A 4), the forward scattering 
sum-rule refers only to wave vectors on the Fermi surface, not to the local version of the Pauli principle. Furthermore, 
the validity of this sum rule has been questioned, pi] The effective mass at this level of approximation is the bare 
one, as in a transitionally invariant local Fermi liquid. p8| Recall however that our microscopic calculations are not 
phenomenological: they explicitly give a value for the Landau parameters. Also, our results extend well beyond the 
q = quantities usually considered in Fermi liquid theory. 

The quasi-particle weight Z calculated with EJj can differ substantially from the initial one. This means that if we 
were to calculate the susceptibility with the corresponding frequency and momentum dependent irreducible vertices 
r*^^) there would be sizeable compensation between vertices and self-energy because our calculations with S^ (Z = 1) 
and constant renormalized vertices already gave excellent agreement with Monte Carlo simulations. 

Finally, consider the high-frequency asymptotics. Since we use bare propagators, the high-frequency asymptotics 
comes in at the appropriate frequency scale, namely ifc„ ~ W and the Hubbard bands do exist in our th eory . However, 
the coefficient of proportionality in front of the asymptotic form 1/ikn is incorrect. Using Eqs.(E9),(Al7),(Aia) we 
can write the high-frequency asymptotics in the following form 



lim 

ikn — *o 



Eo- (k,z/c„) = Un^c 



U 



Us. 



Uch 



^ {nl,) - Uchnl, + (^ 



Usp 



Uch 



\"T"i; 



(ElO) 



This form is useful to understand what is necessary to obtain the quantitatively correct high-frequency behavior. 
Indeed, one would recover the exact result Eq.(p8|), if one were to take into account that: i) the irreducible vertices 
become equal to the bare one U at high-frequencies; ii) the local Pauli principle (fitg ) — n^^ is satisfied. Contrary 
to most other approaches, our theory does satisfy the local Pauli principle Eq.(A15) exactly. However, since our 
irreducible vertices are constant and tuned to describe the low energy physics, we violate the first of the above 
requirements. It is thus clear that for a correct quantitative description of both the low energy physics and the 
Hubbard bands one needs to work with frequency-dependent irreducible vertices. 
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